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Abstract 

We elucidate aspects of the one-loop anomalous dimension of so(6)-singlet multi-trace oper¬ 
ators in A/" = 4 SU{Nc) SYM at finite W . First, we study how 1/iVc corrections lift the large W 
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coefficient and they mix only among the operators with the same number of traces at leading 
order. Second, we study the lowest one-loop dimension of operators of length equal to 2Nc . 
The dimension of such operators becomes more negative as W increases, which will eventually 
diverge in the double scaling limit. Third, we examine the structure of level-crossing at finite 
Nc in view of unitarity. Finally we find out a correspondence between the large N^. zero modes 
and completely symmetric polynomials of Mandelstam variables. 
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1 Introduction 

Large N^ gauge theories have been intensively studied for decades since the work of ’t Hooft, 
who discovered that mesons are described by a string in the planar limit jT]. Maldacena argued 
that a d-dimensional gauge theory can be described by a (d + l)-dimensional gravity theory, 
which is now called the AdS/CFT correspondence |2]. 

The primary example of the AdS/CFT correspondence is the one between maximally super- 
symmetric Yang-Mills theory in four dimensions {M = 4 SYM) and superstring on AdSs x S® 
in the planar limit. Integrability is a powerful technique in this setup, which makes it possible 
to compute conformal dimensions of gauge-invariant operators in A/" = 4 SYM and the energy 
of AdSs X superstring states at any value of the’t Hooft coupling (see [3] for a review). In 
contrast, the non-planar problems are notoriously complicated and less well-understood. The 
entire problem seems daunting, but we shall present a tractable corner of it. The purpose of 
this paper is to look for simple structures in the non-planar problem and to provide a new point 
of view on the large Nc problems. 

Our motivation is explained in three ways. The first one is to better understand the 
AdS / CFT correspondence at finite A7 . 

The second one is to study operators with negative anomalous dimensions. Since all planar 
anomalous dimensions are non-negative, negative anomalous dimensions are inevitably a non- 
planar effect. In addition, 1/Nc corrections to the dimension of double-trace operators are 
related to four-point functions. The four-point functions on the gravity side are related to the 
phase shift of a high-energy scattering ISEIE], which must be positive in order to preserve the 
asymptotic causality of spacetime |7]. Thus, the AdS/CFT and causality predict that certain 
double-trace operators should have negative anomalous dimensions. We want to observe similar 
phenomena in A/” = 4 SYM, though the operators of concern are different from [7|. 

The third one is to inspect the large N,, behavior of determinant-like operators with all 
non-planar corrections taken into consideration. It was found in that, if we take the naive 
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’t Hooft limit and apply integrability methods, then the dimensions of the donble-determinant 
operators, which correspond to the energies of a pair of open string tachyons between a D-brane 
and anti-D-brane, becomes either divergent or complex at two-loops. Since any operators of 

= 4 SYM have pertnrbatively real and finite anomalons dimensions at hnite Nc , snch 
pathological behavior shonld be remedied with non-planar effects. 

Given this sitnation, it is worthwhile to look for an alternative to integrability to stndy the 
non-planar problem. A promising approach is the gronp theoretical method initiated by [TU]| . 
This method is based on the constrnction of a basis of local operators that diagonalizes the free 
two-point fnnctions. The work nm dealt with the 1/2 BPS sector, and later it was extended 
to more general sectors. The basis of mm is snitable for describing open string excitations 
on giant gravitons, and the bases in [131 HIES] were designed by Braner algebras. The flavor 
symmetry is manifested in the bases of [T^ m- One common featnre of the diagonal bases 
is that they are labeled by a set of Yonng diagrams. It was also observed that the one-loop 
mixing is highly constrained on these bases — two operators can mix if their Yonng diagrams 
are related by moving a single box of the Yonng diagram dSl ESI 1231 Ej. This property makes 
the problem tractable, and in certain sitnations it is possible to obtain the spectrnm of the 
dilatation operator explicitly |22l ESI El ESI EHl ETj. Althongh we have realized these bases 
are nsefnl tools for the non-planar mixing problem, most resnlts have been limited to the 5u(2) 
sector. One motivation of this work is to explore the non-planar mixing problem of the so(6) 
sectorfl 

The non-planar dilatation operator in the one-loop so(6) sector was written down in [281129] . 

Sone-ioop = ^ : (-^tr[<h^,<h„][cE)-,<|)T - itr[<h^,<l>’^][<h^,<l>’^]^ :, (1.1) 

and the general non-planar spectrnm has been stndied in many ways [301131] . In the planar 
limit, the Bethe Ansatz Eqnations give ns the spectrnm of all single-trace operators in the 
so(6) sector |32|- Natnrally, one hopes to know the spectrnm of all mnlti-trace operators at 
arbitrary . This problem has not been stndied thoronghly. One of the main diflicnlties 
is that the nnmber of operators involved in the mixing grows factorially with respect to the 
operator length L. For example, there are 469 scalar so(6) singlet mnlti-trace operators at 
L = 10, and 4477 operators at L = 12. By brnte-force compntation, we managed to compnte 
explicitly the matrix elements of the non-planar mixing and obtained their eigenvalnes np to 
L = 10. Varions interesting properties of onr resnlts are explained below. 

At Nc = oo, the anomalons dimensions are highly degenerate, and this is why the planar 
mixing problem is integrable. Most of the large Nc degeneracy are lifted by 1/Nc corrections. 
We focns on the degenerate eigenstate having the zero anomalons dimension at large Nc , which 

^The su(2) sector consists of gauge-invariant local operators made out of two holomorphic scalars, which are 
closed under the operator mixing to all orders of perturbation theory in W = 4 SYM. The so(6) sector is made 
out of six real scalars and closed at one-loop. 
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we call the operator submixing problem. The submixing problem at 0{1/Nl) can be solved by 
diagonalizing the Hamiltonian 

-^S°m ~ limPo ~ T >1 ( 2)0 + e) 2)i] , 2)one-loop = S)o + ) (1-2) 

£^•0 iVc 

where Po is the projector to the space of the large zero modes. We observed remarkable 
patterns among the solutions up to L = 10 , and conjecture that these patterns continue for 
any L. In particular, each eigenstate is a sum of multi-trace operators with the same number 
of traces, and all eigenvalues are non-positive: 

[ ^sm , # (traces) ] = 0 , 72 < 0 , ^ ( 1 . 3 ) 

where 7 = X]n=o In is the one-loop anomalous dimension. 

Besides the dilatation operator, higher-point correlation functions offer an alternative route 
to compute the 1 /W corrections to the operator dimension or string energy. For example, 
by studying the OPE expansion of the four-point functions of single-trace operators, one can 
compute the 1 /W correction to the dimension of intermediate states, which are double-trace 
operators at large ■ This problem has been explored in the literature on the AdSs x side 
j38| IHU |35l [36], [37], the SYM side jSH] EH] [10] E], and via the conformal bootstrap jl 2 ] H3] 
Si 05] 06). Following this line of development, we consider four-point functions of multi-trace 
operators in Section and argue that they constrain the eigensystem of submixing matrices. 

In Section we discuss the following aspects of the hnite W spectrum: 

• The one-loop dimension of determinant-like operators generally receives non-planar correc¬ 
tions even in the large W limit, because their operator length is of order W. The proper way 
to take the large W limit is to study their dimension at a hnite W ~ -h and extrapolate the 
results to W 3> 1. This result can be different from the naive limit, where we hrst take the 
large W limit and identify a double-determinant-like operator as a state of an integrable open 
spin chain with boundaries. We are interested in the operators with length L = 2Nc , with the 
hope that non-planar corrections rescue the pathological behavior of the double-determinant 
operator found in [S]. Indeed, we hnd an operator with L = 2Nc whose one-loop dimension is 
zero at any , which may be regarded as a non-planar completion of the double-determinant 
operator. In addition, we also hnd an operator with a lower dimension by studying dilatation 
eigenstates at W = T/2 up to L = 10. The lowest-energy state has a negative one-loop dimen¬ 
sion which decreases as iV). increases. This anomalous dimension will eventually diverge in the 
large limit 

A - L -)■ -cx), gyu -t 0, W = ^ -t 00 , A = Ws'ym ^ . (1.4) 

^Of course, the tree-level dimension of the double-determinant is 2Nc, which diverges in the large W limit. 
Here we argue that the anomalous dimension in the so (6) sector diverges. This situation is different from the 
one in the su(2) sector, where the determinant operator det ^ dual to giant graviton is BPS [17], and the 
determinant-like operators have finite anomalous dimensions in the large W limit and are computed by the 
energy of an integrable open spin chain |48l H5] . 
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• We investigate level-crossing: whether the adjacent one-loop dimensions or energy levels cross 
at finite W- According to the non-crossing rule of von Neumann-Wigner |50| . there should be 
no level-crossing if the Hamiltonian of the system has no extra symmetry. We find that most 
of the energy levels do not cross for Nc G [L, -|-oo] but do collide for Nc < L, which can be 
explained by the non-Hermiticity of the operator mixing matrix and the finite W constraints. 

• The one-loop dimensions exhibit eccentric behaviors for W < L, which makes it difficult to 
keep track of the dilatation eigenstates (or A/" = 4 SYM operators) from large W to small W. 
In other words, when we regard the one-loop dimension as an analytic curve on the ( 7 , W) 
plane, we are not able to tell which curve corresponds to which SYM operator unambiguously. 
To support this idea, we construct an automorphism which relates high energy states and low 
energy states. 

Finally, Section]^ is devoted to Discussion and Outlook. 


Appendices are organized as follows. Our notation is explained in Appendix]^ The one-loop 
spectrum at finite Nc is summarized in Appendix The details of computations on correlation 
functions are given in Appendix]^ The intriguing relation between the number of so(6) singlet 
large Nc zero modes and the completely symmetric polynomial of Mandelstam variables is 
discussed in Appendix D A concise way of describing multi-trace operators is introduced in 
Appendix where we briefly explain our idea behind Mathematica implementation. This 
paper is accompanied by a Mathematica notebook and a data file which contain the operator 
mixing matrix at L = 10 discussed in Appendix |B.2[ 


2 Operator (sub)mixing problem 


2.1 Synopsis 


We start with a brief review on the operator mixing problem. In A/" = 4 SYM at weak coupling, 
the two-point functions of scalar operators behave a^ 


(Oa(x)Os(O)) 



Sab + T^ab log(|a^l/^) 


| 2 A( 0 ) 

x\ 


+ C>(5'ym)) 


( 2 . 1 ) 


where and T^^b ^(S'ym) come from finite and singular counterterms, respectively 
|4n) . The quantity T^^ = related to the one-loop mixing matrix, or the one-loop 

dilatation operator. By using the eigenstates of T^) we can rewrite the two-point function fl2.ip 
into the familiar form 


{dA{x)dBm ^ 


Sab 


2(aW+aOA ’ 


\x 


= Oa. 


^Note that 0 in logarithmic CFT. 


( 2 . 2 ) 
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Note that F is generally not Hermitian even though S,T are Herniitian j5T[[29]. 

The dilatation operator is self-adjoint with respect to the two-point functions, which can be 
shown by the spacetime translational symmetry as 

(-iT){{C>^(x)OB(0)>-{O^(0)C>B(-a:)>} = 0 ^ {(SOa)Ob) = {Oa(SOb)). 

( 2 ^ 3 ) 

If Oa , Ob are the eigenstate of the dilatation operator D, we obtain 

(Aa — As) Sab = 0 Sab = 0 or Aa = As . (2.4) 

Below we investigate the case A a = As = C>(l/iVc) and Sab 7^ 0, where the spectral degeneracy 
among the large Nc zero modes is lifted by 1/Nc corrections. This problem is what we called 
the operator submixing in Introduction. 

We want to diagonalize the submixing matrix. This problem may look simple because 
the dimension of the submixing matrix, namely the number of the large N^. zero modes, is 
significantly smaller than the number of all operators. Nevertheless, it contains rich information 
about non-planar interactions. 

The submixing equation itself determines the eigenstates at the zeroth order of l/N^. expan¬ 
sion. We focus on the problem at 0{1/N^), because in A/" = 4 SYM, the degeneracy of the large 
Nc zero modes is mostly lifted at this order. Still, some eigenvalues may remain degenerate. 


2.2 Lifting the large Nc degeneracy 


Let us derive submixing equations. We are mostly interested in the submixing among the large 
Nc zero modes, namely the operators with the vanishing anomalous dimension at W = oo. The 


large W zero modes in the scalar so(6) singlets will be given explicitly in Section 2.3 


Let Ho and H, be the space of the large W zero modes and non-zero modes, respectively. 
We split the dilatation operator into the planar and non-planar part as 2)one-ioop = 2)i, 

and regularize the planar part by = e -f 2)o to make it invertible. Take a general linear 
combination of the large W zero modes, and call it 'ipo. 'ipo is regarded as the large W limit 
of a finite W eigenvector ip. Then the non-planar operator mixing equation 2)one-ioop'0 = 7t/’ 
becomes 


OO 


OO 


(I), + 7Va'Di)^= e + 5^iV-7*U, P^ = Y.N;%, (V^o e , 7o = 0). 

\ i=0 J i=0 

Expand this equation in N~^ and apply 2)7^• The results are 

n 

lAn = y 7n-iSr'’/’i + > (" > 0. <^-1=0). 


(2,5) 


( 2 , 6 ) 


j=0 
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We take the limit e —)■ 0 and keep track of the singular terms. Denote the matrix elements 
of 2)7^ and Di by 




ij 




(Oo/. e 'Ho/.), 


(2.7) 


where = A, + e are the non-zero eigenvalues of 2)^. The projector to Ho is written as 


Po = lini e2), \ = Po . (2.8) 

e-s-O 

The projector to H, is P, = 1 — H, • We also use the notation = Po/, tjj. 

The equation (12.61) can be written as 

^ ® 1 ■ (2.9) 

The limiting behavior of the operator depends on the operand; it is singular on T-Lo and 
regular on P,. By evaluating fl2.9p up to second order and using 70 = 0 and = '^o , we 
obtain 


ro = 0, (2.10) 

^* = e“Si V’o - 2)7^2 )ii7o , (2 .11) 

^72 = e“S2V’o - 2)7^2 )i^7i . (2 .12) 

By extracting the components of Ho/, in fl2.1ip and fl2.12l) . one hnds 

= Ho2)7^2)i V’n-i, V’* =-H, 2)7^2)iV'n-i, (n = l,2). (2.13) 


The hrst equation determines the eigenvalue 7^, which may lift the large N^. degeneracy. The 
second equation determines "0* but not 'ip^ . No equations can constrain '0)) because are the 
zero modes of the planar dilatation. They cannot be hxed by perturbation at n-th order. These 
zero modes are not important at 0{N~‘^)^ because the eigenvalue 72 does not depend on 'ipl as 
shown below. 


2.2.1 Matrix elements 

Let us rewrite the above equations in terms of matrix elements. We use small or capital letters 
in place of o or • to label the components of Ho or P,, respectively. An eigenstate is expanded 
as ipn = Un,a Oo,a + Un,A 0,^A ■ Equatious fl2.13p become]^ 

e"^ {mab Un-i,b + rUaB Un-1,B) + 0{e^) 

A7^ {niAb Un-l,b + rUAB Un-1,B) + 0{e) 

'^Note that ip and 7 do not depend on e. 



^ 7r* ^0,a 


'l^n,A 
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Collecting the terms leading in the limit e —?• 0 at n = 1, we obtain, 


0 = {jlSab - rriab) Uo^t , 0 = Mo,a • (2.15) 

Aa 

These are the equations for degenerate perturbation at first order. In Appendix |B.2| , we observe 
that the large Nc degeneracy of zero modes in A/" = 4 SYM is not lifted at first order. Therefore, 
we assume 

7i = mab = 0. (2.16) 


At n = 2, we find 


0 

0 


^aB^Bc 

A_B 


fnABfnBc 

^A^B 


T 72 ^0,c 5 

rriAb 

Uo,c -7- Ui^b — U2,A 


(2.17) 

(2.18) 


The first line is the submixing equations because it lifts most of the large Nc degeneracy at 
0{N~‘^). This equation is insensitive to the ambiguity of 'ipl- The second line determines the 
eigenstates 


2.3 Large Nc zero modes 


We classify the large Nc zero modes in the singlet representation of so(6). 

The BPS states of W = 4 SYM can be labeled by the irreducible representations of 
psu(2,2|4) j52] (see also |53]). The conformal part of the algebra psu(2,2|4) can be neglected 
if we focus on primary operators. Let us denote by the single-trace operator consisting of 
^ so(6) scalars whose flavor indices are traceless and symmetric. General half-BPS multiplets 
belong to the representation [0, L, 0] (L > 2). The primary operators can be written at large 
Nc as 


C^half — 




z=l 


—T, 1/ > 2, 


(2.19) 


J [0,L,0] 


2=1 


Similarly, quarter- and eighth-BPS primary operators can be written at large Nc as 




quarter 






eighth 


, 2=1 

m 


lie' 


2=1 


[h^k^h] 


[h,k,h+2h'] 


li = k + 2h, h > 1, 

2 

ii = k + 2h + 3h', h' > 1. 

2 


It is not easy to write down BPS operators explicitly at finite Nc because one has to solve the 
mixing problem coming from the color structure |38l [5H [55] . Recent attempts at solving the 
zero eigenvalue equation of the dilatation operator are available in j56l (57] 158] |23l l25] . 


®Here [p, q, r] denotes Dynkin labels of su(4)h = so(6). 















One can prove that a large Nc zero mode is always written as a prodnct of At large 

Nc, the anomalous dimension of a multi-trace operator is the sum of the anomalous dimension 
of the constituent single-trace operators. The dimensions of the single-traces are given by the 
planar dilatation, which is equal to an integrable spin-chain Hamiltonian. The planar dilatation 
has the ground state given by the traceless symmetric single-trace operator, which is and 
its eigenvalues are non-negative. Note that the products of BPS operators are non-BPS in 
general. 

We present examples of scalar multi-trace operators in so(6) singlet with the vanishing 
anomalous dimension at large W. For clarity, we introduce the flavor indices 

= tr (<F(q<Fi 2 • • • (2.20) 


where (A *2 ■ ■ - ii) means traceless and completely symmetric. Then, the large W zero modes 
with length L = 4, 6, 8 are given by 


{Cij f'p}) 

{CijkCijk, CijCjkCki}, ( 2 - 21 ) 

{^ijkl ^ijkl ) ^ijkl ^ij ^kl y ^ijk ^ijl ^kl y ^ij ^ji ^kl ^Ik y ^ij ^jk ^kl y 

The number of the large W zero modes with length L, denoted by Zl , grows as in Table 


L 

2 

4 

6 

8 

10 

12 

Zl 

0 

1 

2 

5 

11 

34 

dim^L 

1 

4 

15 

71 

469 

4477 


Table 1; dimT^L counts the numbers of all scalar multi-trace operators in so(6) singlet with 
length L. counts the operators in whose anomalous dimension vanishes at large W. 


Finding an explicit formula for Z^ is a difficult problem of combinatorics. Interestingly, this 
series coincides with the asymptotic number of completely symmetric polynomials of Mandel¬ 
stam variables at degree L/2 subject to the massless momentum conservation [3^ [UU] , which 
will be further explained in AppendixThis coincidence may break down at L > 14 due to the 
finite-A^/ constraints, where A^j = 6inA/' = 4 SYM. An example of the finite-Y/ constraints 
is the following anti-symmetrization identity: 

0 = sign((j) ^* 2 V( 2 ) ^*3V(3) ^*4V(4) ^*5V(5) ^*6V(6) ^*7V(7) (*A: = 1, 2,..., 6), (2.22) 

o-giSy 

which reduces the number of independent large zero modes by one. 


2.4 Observations on submixing 


We present here interesting structure found in the submixing. The data in Appendix B.2 


IS 


summarized in Figure [T| which shows the eigenvalues of the submixing Hamiltonian for 
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L/2 

2 

3 

4 

5 

-72 

20 

60.2459 

101.276 

132.315 


Table 2: The coefficient 72 of double-trace operators for various L. 


K2 



Figure 1; Plot of 72 for the large zero modes against the operator length L at iVj = 6 . 

L = 4,6,8,10. The one-loop anomalous dimension is given by 7 = 72 /fV^ + 0{N~^). The 
legend “m-tuple trace" means that the operator consists of a mixture of m-trace operators at 
large . At higher orders, they mix with everything else. Also, the numerical values of 72 for 
the operators starting from double-traces are shown in Table 

Let us rephrase our findings. First, most of the large Nc degeneracy of zero modes is lifted 
at second order. This statement is not trivial because the large Nc degeneracy of non-zero 
modes is sometimes lifted at the first order in l/Nc- 

Second, all the second-order eigenvalues are non-positive, 72 < 0. This is again non-trivial 
despite the minus sing in fl3.30p . because the matrix elements are not symmetric, maA 7 ^ ^Aa • 
Note that the finite Nc anomalous dimensions of the large Nc zero modes can be positive. There 
exist operators with 70 = 7i = ■ ■ ■ = 7n-i = 0 and 7n > 0 for some n >3. 

Third, when the large Nc degeneracy of zero modes is lifted at 0{1/N^), all eigenstates 
have a definite number of traces. Thus, the submixing Hamiltonian should commute with the 
operator which counts the number of traces. In contrast, the operators with different number 
of traces may submix for the large Nc non-zero modes. 

Fourth, if we define the density of submixing matrix elements by applying the Laplace 
transform to 2 )g ^ a^ 

Kmi't)ac = - y] UIbc , [ dt hac{t) = ^ (2.23) 

^ Jo 

®There can be other integral representations of 2 )q 


10 












the submixing densities at different t satisfy 

0 = -Pmin )fec ~ ^sm(^ )abhs^{t)bc^ 

~ abhs^{t )bc — h^^(t )abh^^{t)b(^ Pmin ~ 0) (2.24) 

where is the projector to P^in, which is the subspace of Po spanned by the products of 
length-two operators. This “projective commutation” relation is equivalent to 

0 = Pmmi^aB (Afi) mscmcD (Afl) I^De “ (A^) mBcrricD (A^) niDe^ 

= (niaB (As) rriBcmcD (A^) moe - rUaB (A^) niBcmcD (A^) rriDe'^Pmm, (2-25) 

for n,m>0. These relations imply that the submixing problem on is integrable. 

It will be interesting to prove these findings and see how general these relations hold true, 
which will help to solve general submixing problem. The projection by P^j^ may be removed 
by finding a better definition of the submixing density, or equivalently a better definition of 
mutually commuting charges fl2.25p . 

2.5 How to solve non-planar operator mixing 

Let us pause to explain how we obtained the aforementioned results on the operator submixing. 
It is not a priori clear at which order the large N^. spectral degeneracy is lifted, so we need to 
study the general non-planar operator mixing problem. Since the full non-planar problem is 
too complicated to work by hand, we used computer programs extensively to obtain concrete 
results. More specifically, we wrote Mathematica codes and proceed as follows, 

• Generate a complete set of multi-trace operators with a given length, using the notation 
in Appendix [E} 

• Compute the action of one-loop dilatation operator on the complete set of multi-traces. 

• Extract the matrix elements, as analytic functions of Nc and A^/ 

Our computation can be done symbolically with no numerical errors until this stage. See the 
attached Mathematica code for details. 

Having obtained explicit non-planar mixing matrices, we will study their eigensystems later 
in this paper. We will analyze 

• Operator submixing problem as 1/W perturbation 

• Eigenvalues as analytic functions of Nc and Nf 

• Finite W constraints on eigenvectors 

There is no need to invent new techniques to solve these standard problems of linear algebra. In 
order to manage large matrices (of size < 5000) efficiently, one should compute their eigenvalues 

= 4 SYM corresponds to Nj = 6 . We introduced a parameter Nf as discussed in Appendix IaI 
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numerically at MachinePrecision. It will take an extremely long time to find all eigenvalues 
analytically if the matrix size is bigger than ~ 30. 

It is not easy to inspect carefully the Mathematica codes, so we need to check that no 
mistakes were made during the computation, which can be done in two ways. First, our mixing 
matrix eigenvalues at L = 4, 6 agreed analytically with the literature ISH 162 ] as shown in 
Appendix B.2 Second, all eigenvalues of the mixing matrix are real. This is a non-trivial 


check because we computed the mixing matrix elements using the basis in which the matrix is 
non-hermitian. Further discussion on the hnite Nc spectrum will be given in Section 


3 Correlation functions 

We will study the submixing problem by evaluating two-point functions. First we evaluate 72 
in terms of correlation functions of 'ipi, and derive the submixing Hamiltonian. Then we discuss 
the sign and the L-dependence of 72 . 

3.1 Preliminary 

In this subsection, we summarize preliminary facts about correlation functions. 

We will focus on color and flavor combinatorial factors of correlation functions, and we 
drop the space-time dependence, which is trivial to recover from conformal invariance. In our 
analysis, the difference between the U (W) gauge group and the SU (W) gauge group does 
not matter because we mainly discuss the leading large W contribution. We use the U{Nc) 
Wick-contraction 


It can be extended to the case of multi-helds as 


(3.1) 



■ (41,. 

... Ai 

/mi K^OLJrriL 



o-£Sl 


Ja(l) 

"C. 


= ^aiba(l) ■ 

^a.Lb,r(L) 

>J ■ 


(3.2) 




where the normal ordering defines an operator without self-contractions. Because of the normal 
ordering. Wick-contractions apply only between and . 

Next gauge invariant operators are considered. Any multi-trace gauge invariant operator 
built from L matrices can be characterized by an element of the symmetric group Sl as 


trar4f)^(r);(4r)i 


J (^®L\I 


*r{l) 






(3.3) 
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Focusing only on the color structure by leaving out the flavor indices, this description has the 
symmetry 




(3^4) 


where g is any element in Sl. This implies that the multi-trace structure of gauge invariant 
operators is classified by the conjugacy classes of Sl, which are the partitions of L. It is then 
straightforward to obtain two-point functions for gauge invariant operators DU, 

CT&Sl 

= • ■ ■ <5aiC^(i)trL('rin T2(t) 

CT&Sl 

= (3.5) 

CT&Sl 

We introduced a shorthand notation 




■•■5a 


= (a|cr|c), 


(3.6) 


and C{a) counts the number of cycles in the permutation a, e.g. C'(l) = L, C'((12)) = L — 1. 

Finding the L-dependence of the correlator fl3.5p is a complicated combinatorial problem 
with the color and flavor structures correlated. In Section 3.3, we will concretely compute 
the two-point functions of the double-trace large W zero mode, as an example to see how the 
correlators depend on L. 

In the literature nanziiii], bases of gauge invariant operators in the so (6) scalar sector have 
been constructed, and they diagonalize the two-point functions including all 1/W corrections. 
The readers can also refer to Uni ESI [HI E3], although we do not need those technologies in 
what follows. 

When Ti and T 2 are in the same conjugacy class (i.e. they have the same trace structure), 
the two-point functions behave as ~ N^. Its coefficient is a symmetry factor determined by 
the number of permutations p satisfying Tip~^T 2 P = ij^ 


{: ) :: tri(T24f) :> = 5^(a|p|c> + O(Nt'). (3.8) 

P 


® The symmetry factor of the two-point functions can be associated with group theory. Suppose t is in the 
conjugacy class c = [ci, • • • , cl], where is the number of cycles of length i, i.e. L = ici. The number of 
permutations p satisfying p~^Tp = t is counted as follows. For a cycle with length i, there are exactly i cyclic 
permutations that do not change the cycle. Hence, cyclic permutations can leave the cycles. Also, p can 
permute Ci cycles of the same length. Then we have the formula 


SyrriLic) = dim 


0(5,®Zf-) 


L 




(3.7) 


For T in the conjugacy class cl = 1, we have SymL{c) = L. For r in the conjugacy class cl /2 = 2, we have 
SyrriLic) = I? 12. 
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On the other hand, if Ti and T 2 are not in the same conjugacy class, (i.e. they have different 
trace structure), the two-point functions behave as with a < L. In other words, for two 
normalized operators (pi, ip 2 we have 


{1P1IP2) ~ 0 ( 1 ), 

if they have the same trace structure, and 

{tpitp 2 ) ~ (a > 1). 


(3.9) 


(3.10) 


if they have different trace structure. 

We can also deduce the large W behavior of matrix elements of the dilatation operator. 
The planar dilatation operator does not change the number of traces while the non-planar 
dilatation operator does change the number of traces by one. In other words, the planar 
dilatation operator does not change trace structure of operators while the non-planar operator 
does. As a consequence we hnd 

0(1), i^^~o(iv;‘), (3. 

and so on, where 0 is a linear combination of operators with the same trace structure. 

We expand a dilatation eigenstate as 

1 1 

= + (3.12) 

n=0 




where each held 'i/'i is allowed to be a linear combination of multi-traces. We will assume that 
the two-point function between ijji and 'ipi+i is subleading 


(VaA+i) 


~ 0(iV-^). 


(3.13) 


For example, if "00 is a double-trace operator, 'ipi can be a linear combination of single-traces 
and triple-traces. This assumption is essential in the subsequent argument, and it is desirable 
to have a proof. 

Note that correlation functions of '0 like (13.Sh have a meaningful 1/W expansion only if 
Nc S> L. In general, the numerical factor appearing in each order of the expansion is a function 
of L 01 Li, where Lj is a partition of L. Hence the expansion parameter in fl3.13p and fl3.1ip 
is f{Li)/Nc. In order for this expansion to be a perturbative expansion, we need to assume 
Nc 3> L, which we do in what follows. Then, two-point functions of 'ijj reduce to those of V’o in 
the planar limit. 
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3.2 Correlator expressions of 7 

3.2.1 1/A^c perturbation revisited 

The dilatation operator can be divided into the planar and non-planar parts, 

2)one-loop = 2)o + 77 (3-14) 

We will study the case V’o is given by a linear combination of the large W zero modes, which 
is characterized by 

DoV’o = 0. (3.15) 

Substituting the expansions into the eigenvalue equation 

2) one-loop^ = 7'0, (3-16) 

we have 

^ (2)o d'lp + Di -00) + '5V' = 7 (^^0 + ^ ^V'^ • (3.17) 

Taking the two-point functions of this equation with -00 leads to 

^ (V’o (2)o<5V' + SiV'o)) + ^ (V'o (Si'^V')) = 7 (^(V’oV'o) + ^ (V’o<5V')^ • (3.18) 

From the fact that the dilatation operator is self-adjoint with respect to the two-point function, 
we have 


(V’o i'OoSi’)) = ((SoVo) VV) = 


(3.19) 


Because (J'V = V'l + V’ 2 /lVc + ■ ■ ■, we obtain 


Ncj'lpO^l'lpo) + (V^pDl-Vl) 

mo) ^ ^ 


(3.20) 


Assuming that V^o and Hiv’d have different trace structure, we see that the expansion of 7 starts 
from the order Substituting fl3.20p into (13.17p . we find 

HoV’i + HIV’D = ^, (3.21) 


where 7 is a linear combination of multi-traces that are different from ipi. Considering the 
two-point functions of (13.1711 and V’l, we obtain 


(V’iV’d) = ^^((V’w) + (V’lSiV’i)) - ^ (V’lV’i) = 0{N^ ^)\/(V’dV’d) (V’lV’i) , (3.22) 

which is consistent with our assumption (13.1311 . 

®This situation is parallel to our argument around (I2.16|) . 
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The equation (3.21) can be used to get the following expression, 


7 = 


1 Nc (-00 - ('01 2)o'0i) 


+ o(iv-^), 


(3.23) 


(00 0 o ) 

or equivalently by expanding 7 = J 2 n=o ^c^ln, 

„ „ A^c (0o2)i0o) - (0i2)o0i) 

7o = 0, 7 i = 0, 72 = - - 

(00 00 ) 

This expression of 72 will be frequently used later. Because 72 is the leading term of the 1/Nc 
expansion, we can use the planar limit to compute the correlation functions. 


(3.24) 


We will next give another form of 72. Following the methods of Section 2^, it is not difficult 
to rewrite the above expression of 72 as 


(00 (iVcSi-2)12)712)1) 0o) 

72 = lim- J-—— - 

6^0 (00 0 o) 


This simplihes further as 


72 = - lim 


(T2),T) 


(00 0o) 


T = 2), i2)i0o - -^00, 


(3.25) 


(3.26) 


which can be checked by using the self-adjoint property of 2), and (3.15), (3.21). 


3.2.2 Submixing Hamiltonian 

Let us derive the submixing Hamiltonian. We use a Greek letter to label different large W zero 
modes as 

2)one-loop0^“^ = 7^“^ = 0. (3.27) 

For simplicity, we assume that the eigenvalues are non-degenerate, 7*-“) 0 7^^) for a ^ (3. By 
generalizing the above argument, we obtain 

72“^ 0 S“^) = Nc (^ 0 ^^^ 2 )i 0 ^“^^ - ^ 0 f ^ 2 )o 0 i"^^ 

= lim (4^\Nc^i - 2)i2)7i2)i)0i“^\ . (3.28) 

The left-hand side is proportional to 6 ^ 13 , so the right-hand side is. Since {00°"^} spans the 
entire Hilbert space of the large W zero modes, we can think of fl3.28p as an operator identity, 

4 “’= 7777 . = ( 3 . 29 ) 

which is fl3.3ip . 

In Section we will rewrite the submixing equation in the operator form, 

72 Mo, a Oo,a = “ Uq^c Oo,c < -1 72 00 = 0o , (3.30) 

= Pc i7sm = Po (w 2)i - 2)i 2)71 2)i). (3.31) 
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We call i/sm submixing Hamiltonian. The projector Po is needed because H^ui maps "Ho to the 
whole space of operators, He 
down matrix elements as in Section 


H,. Note that the first term in fl3.3ip vanishes if we write 
owing to (I2.16p . The (dis)appearance of the first term 


2.2 


originates from the fact that the mixing matrix is not Hermitian whereas the dilatation operator 
is self-adjoint with respect to the two-point functions. 

We can also rewrite the submixing density into an operator form as 




'N 

hsm{t) = ( ^ 2)e 


2)l 






2)i 


(3.32) 


The integral over [0, cxd) converges if we consider Pohsmit) = thanks to lime_s.o Po = 0 

and fl2.16p . Then we can safely take the limit e —)■ 0. 

For the double-trace operator CgC ^, the action of the submixing Hamiltonian can be schemat¬ 
ically represented as 


H° 


Nr, 



(3.33) 


where each circle represents the single-trace operator Cg. The Si, represented by a gray circle, 
is an interaction vertex and is a propagator. The intermediate state in the second term 
must be a single-trace. 

In the literature, the energy spectrum of string in pp-wave SFT was compared with the 
dimension of BMN-type operators in A/” = 4 SYM. It was argued that we should include four- 
point contact terms in the pp-wave SFT for a better agreement and to cancel a divergent term 
from a cubic coupling 


|. We can interpret the first term in (I3.33p as the four-point contact 
term, and the second term as the cubic term squared, in consistent with the expectation from 
the pp-wave SFT. 

Interestingly, we can remove the four-point contact term if we renormalize the wave-function 
as T' = —2)^4/ = ^'00 — 2)7^2)i'0o- From (I3.26j) it follows that which 

contains only the propagator 2*7^. 


3.2.3 On the sign of 72 

Here is another comment on the sign of 72. The eigenvalues of 2)o are non-negative, but it does 
not mean 72 is non-positive. Let V^,Vq be an eigenbasis of 2),2)o, respectively. The state T 
can be expanded as 


T = F/ W = F/ F'jHo 


(3.34) 
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Note that U is not an orthogonal transformation, becanse {V^^} is not an orthonormal basis at 
hnite ■ It follows that 


(^ 2 ),^) = {{FkU^ jV^^) 

= FkU^l{V,^1^V,^)FjU^j 

= FK{u\^^{U^fK)Fi. (3.35) 

Since 70 and (t/7of/^) have different spectra at hnite Nc , (ft' De ft') can be positive or negative. 


3.3 Two characteristic classes of the large Nc zero modes 


As we mentioned in Section 2.4, the snbmixing pattern can be ctassihed by the number of traces 
of the targe Nc zero modes. We will focus on the two characteristic cases in the classihcation: 
the number of traces being minimal or maximal. The hrst case is the double-trace operator 
Cavac/ 2 ^o.v- 0 ‘Li 2 ’ ^hich is the unique operator of length L with two traces. The second case 
is a set of operators having the largest number of traces, namely L/2-tuple length-2 operators 
such as CijCjkCki ■ ■ ■ Cmi ■ The number of such operators is equal to the number of partitions 
of L/2 excluding the partitions that contain 1. Let us introduce a shorthand notation = 
Caia2Ca2a3Ca3a4''' Capai- There are two operators for L = 10; and and there are four 

operators for L = 12; and For p > 7, the constraint (2.22) has 

to be taken into account. 


Below we want to see how the L-dependence arises from correlation functions, by taking 
the double-trace operator as the simplest example. In general, the two-point functions of other 
large W zero modes result in complicated combinatorial problems. 

We write the singlet double-trace operator as 

XL = : tr (<F(„, • ■ ■ tr (<!>(„, • ■ ■ ^ac/ 2 )) ■ 

= {d\P[L/2 ]\b) : tTLi{aL/2 ® «L/2) ® /') :, (3.36) 

where 0^/2 represents the conjugacy class of cycle length L/2 in Sl/ 2 . We introduced the 
projection operator P[l/ 2 ] associated with the rank-L/2 symmetric traceless representation of 
so(6). 

From the formula ( |3.5[ ), the two-point functions of xl can be computed as 

(XlXl) = '^{d\PL/2\b){c\PL/2\d){d®b\p\c<» (m/2®m/2)p)_ 

a,&,c,d 

To simplify it, let us use the diagrammatic computation presented in |M1113] and Figure]^ We 
also introduce a trace over where Vf is a 6-dimensional vector space for the flavor index. 
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Figure 2 : The flavor structure of (xlXl)- Two lines at the top and the bottom of the figure are 
identified. The contraction operator ]^. Cu connects a with b, and c with d. 


It follows thatE3 


(XlXl) = ® l)p(^[L/2] ® l)CuC22 ■ • • CL/2L/2) ivf 

P&Sl 

= 2 tr{ [{P[L/2] ® l)(pi ® P2 )(T’[l/ 2] ® l)C'iiC'22 ■ ■ ■ c 112112) 

Pl,P2SSi/2 

X Qrp^L-l\ 


= 2Dimi/2 N^ + 0{Nt^) 

= y iVi'Dimi/2 + 0{Nt^), (3.38) 


where tr{ is the trace over The factor 2 at the second line originated in the symmetry 

that exchanges two 0 ^/ 2 ’s of the operator. At the third equality, and p® are permutations 
satisfying pf'^aL/ 2 pf^ ^ = ttL /2 (see footnote 8), and we have used 


f^L((T’[L/2] ® l)(pi <8) P2){.P[L/2] ® l)C'llC'22 ' ' ' ClI2LI2) 

— fr{((T’[i;,/2] ® l)C'iiC'22 ■ • • CLI 2 LI 2 ) 

= ^^i/2iP[L/2]) 

= Dimi/2, (3.39) 


^‘^With the notation V®^ = 0 Ca is the contraction between the i-th space in the former 

and the z-th space in the latter V®^P. 
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^ = P[r/5i. The second equality can be explained by 


where we used piP[L/2] = P[l/2] and • 

Figure |2 without permutation p. Dimj;,/2 = i^(T/2 + l)(L/2 + 2)^(L/2 + 3) is the dimension of 
the totally symmetric traceless representation of so (6). 


3.4 Evaluating 72 


We now evaluate 72 based on the expression (3.24). 

First, we see that the dilatation operator Sone-ioop HI-1)) annihilates the single-trace operator 
Cai,---,as = tr (*^*(01 ® ® *ha7). The first term of the dilatation operator annihilates Caj^,...,as 

because the indices are symmetrized, and the second term annihilates because the indices are 
traceless. Then, the dilatation operator acts on the product CgCj; = Ca^,... ,bt as 

2)one-loop(CaCg) ~ (2? one-loopCa) Cj + Ca (2?one-loop Cg) + {dCs){dC^) = ((9C^) (c?Cg) . (3.40) 

The two derivatives in the dilatation operator do not act on the same single-trace. 

With the above facts, we will turn to the correlator expression of 72. Suppose that the 
dilatation operator acts locally on the factor 


^ — Cai,--- ,c 


(3.41) 


contained in 'ipQ. The effect of the non-planar dilatation operator is to make a single-trace from 


the double-trace, which we find from the formula (A.7)-(A.14), reducing the number of traces 


of '00 by one. The total action of the dilatation is the sum of local actions on each pair of 
single-traces. 

From the intuition of AdS/CFT, we may regard fl3.41l) as a two-string state. A single-string 


state is also created by the action of the dilatation. In this sense, the first term of 72 in (3.24) 


is considered to be the overlap between the two-string and single-string states, and the second 
term of 72 can be interpreted to be the planar energy of the single-string state. 


We next discuss the sign of 72 in (3.24). Recall that the numerical study for L < 10 shows 


that 72 is always non-positive for -00 given by a linear combination of the large W zero modes. 
It is easy to see that the second term of 72 is non-positive because the eigenvalues of the 
planar dilatation operator are non-negative. The sign of the first term is non-trivial because 
the non-planar dilatation operator contains several terms with different signs. 

In order to get more insight on 72, we will analyze the correlation functions exploiting a 
large Nf limit, which is defined by letting the flavor indices run over 1, ■ ■ ■ ,Nf and by taking 
Nf S> L. (In Appendix we will realize that l/Nj appears with a factor L.) This limit 
describes the case Wj = 6 in a good numerical accuracy for small L as will be shown in Figure 
14 Thus, we expect that the anomalous dimension as given by the ratio of correlators is not 
very sensitive to the value of Nf/L. 

Concrete computations are doable at large Nf, as will be given in Appendix]^ In Appendix 


C.l we will show that the hrst term of 72 in (3.24), (0o2)i0o); is non-positive. In Appendix 
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C.2 we will compute 72 for the two classes of the large N^. zero modes presented in the last 
subsection. We write 72^^ for the double trace operator and 72^'^^^ for the L/2-trace operators. 

In general, the first term of 72 is of order Nf while the second term is of order Nj. Only 
the second term, —{'ipiDo'ipi)/{'ipo'ipo), is important at large Nf . For the double-trace operator 


(3.36), our result (C.21) reads 


( 2 ) 

72 


a 


-N}L^ + 0(N,), 


(3.42) 


where a is an iV/-independent constant that may depend on L. In view of the AdS/CFT 
correspondence, the double-trace operator would correspond to a non-supersymmetric two- 
string state, where each string has the angular momentum ±L/2 with the opposite sign. Then 
72 ^ measures the energy difference between the two-string state and that of the bound state. 
As L becomes large, one needs the more energy to break the bound state into the oppositely- 
rotating two-string state. This interpretation is consistent with the numerical results suggesting 
that aL‘^ is an increasing function of L. 


For L/2-trace operators, we have the following result (C.29), 


(T) 

72 


+ ^LsA)oi{T-,s) + 0{Nf), 


(3.43) 


where T specifies an eigenstate among various L/2-trace operators, and a(T;s) is another Nf- 
independent constant. More details are presented in Appendix |C.2[ These operators should 


correspond to bound states of non-supersymmetric L/2-particles rather than strings. 

Finding the exact L dependence of the submixing eigenvalues is a future problem. One 
of the main obstacles is to compute the inverse of the planar dilatation operator The 

eigenstates of 2)o can be readily constructed by using the integrability methods, but one needs 
to solve combinatorial problems to expand the states Di'ipo and "01 in terms of the planar 
eigenstates. Group-theoretical methods will be useful for this purpose. The limit A'"c ^ T 3> 1 
with Nf = 6 is also worth investigation. Our results imply that the expansion parameter, or 
the effective string coupling, in this region is y/aL/Nc for 03.421) and ~ a/ aL/Nc 03.43p 


4 Submixing from four-point correlators 

We discuss the relation between operator submixing problem and four-point functions of multi¬ 
trace operators. 

Consider the four-point function of half-BPS operators, The four-point 

function of A/” = 4 SYM can be decomposed into a product of three-point couplings and 
super conformal blocks jnSj. By expanding the four-point function in Taylor series around 
gyu = 0, we obtain a sum of anomalous dimensions averaged over all possible intermediate 
operators of given quantum numbers, where the weight of average is given by the product of 
three-point couplings. 
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Figure 3: OPE decomposition of a planar four-point function. 


If we normalize the two-point functions as 0{1), the perturbative part of the four-point 
function scales as 0{1/N^) in the large limit |671168] . The three-point coupling scales as 
0(1) if the intermediate operator has the same trace structure as : C^p'>C^p'> and 0{1/Nc) 
otherwise. For the former case, the anomalous dimension of the intermediate operator should 
scale as 0{1/N^). If not, the anomalous dimension scales as 0{1). This structure is illustrated 
in Figure 

The case of p = 2 has been studied in detail in [3D]. From the four-point function 
(n^=i extract the contribution of so(6) singlets which propagate along the 

internal line. There are four so (6) singlets at twist four, corresponding to the L = 4 states 
flB.6|l . If we denote the eigenvalues of the one-loop dilatation operator by { 7 ^^^}, we obtain 




4 


1=1 


1=1 


9ym ^ ^ 
Svr^ 


Ar 2 ’ 


t-r 

1=1 




W’ 

' C 


at Nc ^ 1. The negative sign of originates from the negative-mode of double-trace 

fyP®) • ■ ■ 


General submixing problem beyond the double-trace is related to the four-point functions 
of products of half-BPS operators. A simple example is 


(:C(2)c(2); (a;^) C^^)[X2) C(2)(a;3) : (a;4)> , (4.2) 


whose OPE decomposition contains information about the triple-trace operator : :. 

Although product operators like : : are not dilatation eigenstates at finite N ^., we can 

neglect 1/Nc corrections by taking the large N,, limit in what follows. 

Some multi-trace four-point functions are 0{1/Nl) or less. In such situations, we can 
deduce constraints on the eigensystem of submixing matrices. Consider the multi-trace four- 
point function shown in Figure]^ 


F = (:C(2)c(2);(a;,)c(2)(a;2)C(=')(a;3)C('Ha;4)> ~ C^2,2W 

I 



■^12-^34 

rr.2 ™2 

"''13"''24 


.7.2 7.2 .. 

"''14"''23 \ 

7,2 7,2 y • 
"''13"''24^ 

(4.3) 
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Figure 4: A four-point function of 0{1/Nj^). The right hgure shows its color structure. 


Since F has the color structure of the hve-point function of half-BPS operators, it scales as 
0{1/N^) at large Nc. Thus we hud 

^Io«"„)~EC'p- 2),2,;C3,3,/(7<'')”~C'(1 /'Vc) (foranyn>l). (4.4) 

/ 

which implies 

C{2,2),2,/C'3,3,/7^^^ ~ C>(l/^c) (for any/). (4.5) 

Is there a situation where both 0 (^2, 2 ), 2,1 and are 0(1)? It can happen if an intermediate 
operator is given by a linear combination of the double-trace and triple-trace large N^. zero 
modes, 

= Cl : c(2)c(2)c(2) ; + C2 : C(3)c(3) . _ 

If is an eigenstate of the submixing Hamiltonian, then its dimension must be 0(1/A"^) 

from fl4.5p . Indeed, our computation with so (6) singlets at L = 6 shows that either ci or C2 
vanishes in fl4.6p . and one of the three-point couplings is 0(1/A^c)- Moreover, the triple-trace 
large N^. zero-mode has zero anomalous dimensions at any . 

Certainly, this argument can be generalized beyond so(6) singlets, and to the states with 
L > 6. Typically, a multi-trace four-point function is 0{1/N^) if the total number of traces 
is different between two sides of the OPE decomposition. Then, from (14.51) the large N^. zero- 
modes which propagate the internal line should have a dehnite number of traces to avoid 
0{1) three-point couplings, or their anomalous dimensions should be 0{1/N^). This result is 
consistent with our observation in Section 12.41 

It will be interesting to check this argument by perturbative results. It is conjectured in [TO] 
that the one-loop n-point correlators of the chiral primaries satisfy a generalized factorization 
formula at large N^.. A similar structure may be found for the multi-trace four-point functions 
considered here. 

5 One-loop spectrum at finite and analytic Nc 

In this section, we want to study the hnite N^. physics intrinsically, not by summing 1/A^c 
corrections using the large theories. We will regard an anomalous dimension as a real 
analytic curve on the (7, plane and study the properties of the curves. 
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5.1 Double determinant operators 


Let us start with the hnite problem. In particular, we consider the spectrum of determinant- 
like operators based on our spectral data at arbitrary Nc . Foundations of the following discus¬ 
sion are summarized in Appendix |B| 

When the length of an operator is comparable to N ^., the mixing between single-trace and 
multi-trace operators is no longer negligible even at large W. One cannot a priori guarantee 
that the dimensions of such operators can be computed by integrability methods. In [8], one 
of the authors studied the following double-determinant operators]^ 


o,,5.|z^^y 






^kik2...kpf^ ••• 1 Jiff h ■ ■ ■ In-1^ -'Jiv ’ O- ) 




which should correspond to a pair of open strings ending on giant- and anti-giant-graviton 
branes. The dimension of the YY operator was computed by integrability methods and pertur¬ 
bative A/" = 4 SYM techniques, which precisely agreed for L, L' > 2. However, at hnite values 
of the’t Hooft coupling, the dimensions of these operators exhibit pathological behavior, which 
can be interpreted as the existence of tachyons mu- Also, when L or L' = 1, an unexpected 
divergence is found at two loops. Thus, we should take the large Nc limit more carefully because 
the operator (15.11) is not the exact eigenstate of the one-loop dilatation. 

In the heuristic argument of [8j, it was assumed that the anomalous dimension of the 
YY operator fl5.1ll becomes non-trivial starting from the wrapping order. In particular, its 
anomalous dimension must be zero at one-loop at large Nc . To obtain the desired result by 
solving the non-planar operator mixing, we can look for an operator of length ~ 2Nc whose 
one-loop dimension is zero at any Nc . Since the length of double determinants exceeds Nc , it 
is important to guarantee that such a state survives under the hnite Nc reduction discussed in 


Appendix B.l 


Zero mode at finite Nc . For simplicity, we study double determinant operators without 
insertion of Z^’s instead of the YY operators fl5.ip . 


Oyy = e. 


2122. 




^k\k2...kff ^ 


ji ■ 


yiN-lyiNy^^ 
^Jn-i ^ In ^ h 


y^N-iykN 

In-1 Jn ' 


(5.2) 


which corresponds to a pair of giant and anti-giant graviton branes without open strings. The 
symbol -I-... represents the terms induced by the operator mixing, such as the one given by 
the interchange Y -H- Y or YY -H- XiXj. 

Our goal is to hnd completion of fl5.2|) with zero anomalous dimension at any Nc . Holomorphic- 
antiholomorphic operators are good candidates. They belong to the so-called k = 0 sector of 
na, and their explicit forms are 

Ohol = f ntrr'd (ntrF'd, 1 . (5.3) 

V a=l / Vfe=l Jab 

^^Here Xk = $2fc-i -t i^ 2 k = {W, Y, Z) are complex scalars of W = 4 SYM. 


24 









The one-loop dilatation fll.ip may break the holomorphic-antiholomorphic structure as 


trY^^trY^” ^ |tr tr (Y^^-^XiY^^ ^Xj), tr ^), ...|. (5.4) 

Once broken, the holomorphic-antiholomorphic structure cannot be restored by further dilata¬ 
tion actions. Thus, the operators fl5.3p do not appear on the right-hand side of D-Oj = MjjOj . 
As a result, if an eigenstate of the dilatation operator contains a holomorphic-antiholomorphic 
element, its eigenvalue must be zero: 

-0 ~ C>hoi + • • • ^ 7[0] = 0. (5.5) 

Similarly, the k = 0 operators in the su(2) sector have zero anomalous dimension |57j . 

Consider the finite A), reduction on the holomorphic-antiholomorphic operators. If L is a 
multiple of four, then the following holomorphic-antiholomorphic operator does not become 
null for Nc > 

Ohoi = (tr (F^)) (tr (F')) = 2^/2 + FI2) (F21F12 + F^^^)+ ..., (5.6) 

One can show that the mixing tr (F^)tr (F ) o (tr(FF))^ cannot cancel (15.6p completely. 
Thus the term fl5.6p survives at Nc = 2. And if an eigenstate is non-trivial at Nc = 2, then it 
must not be null for Nc> 2. 


Closed subsector of one-loop dilatation. Let us briefly explain how the dimensions of 
FF operators are related to those of so (6) singlets which we studied earlier. 

Under the dilatation actions, the FF operators like fl5.2p mix with so(6) singlets as in fl5.4p . 
Since the dilatation itself is so(6) singlet, the so(6) singlet operators mix among themselves. 
Let us define the characteristic polynomial of the mixing matrices for two subsectors. 


^vf( 7) = det {Mij - 7 6 ij) 


YY + singlet 


^singlet (7) = det (M/j - 7 5 lj) 


singlet 


(5.7) 


Since the mixing matrix for FF contains the mixing matrix for singlets as a subset, the poly¬ 
nomial ^yy(7) is divisible by ^singlet (7)- 

Our calculation at L = 4, 6 reveals that the most negative eigenvalue belongs to ^singlet (7), 
although ^yy(7)/^singiet(7) coutaius some negative modes. Assuming that this is a generic 
pattern, we will study the lowest eigenvalue among so(6) singlets in detail. 


The lowest eigenvalue at finite Nc . At large Nc , the state with the smallest 72 shown 
in Figure which is the double-trace operator for L < 10, has the lowest eigenvalue among 
all so (6) singlets at a fixed L. Moreover, we observed that this lowest eigenstate does not 
become null for Nc > L/2. The numerical values of the lowest eigenvalue in the so (6) singlets 
dX Nc = L/2 are shown in Tableand Figure]^ 
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L 

4 

6 

8 

10 

7 

-3 

-3.22476 

-3.62944 

-4.01516 


Table 3; The smallest one-loop anomalous dimension Ai = (^)7atiV, = L/2. 


K 



Figure 5: Plot of the lowest anomalous dimension at = Lj^ in Table The quadratic £t is 
given by 7 = -2.70242 - 0.031668L - O.OIOOOL^. 


The results show that I 72 I, the coefficient of the one-loop anomalous dimension, increases 
faster than linearly as increases. If we extrapolate this behavior at large , the lowest- 
energy state will have the dimension 

A = 2iV, - (ao + «iiVc + Oi^Nl) -7 ..., (5.8) 

neglecting a* for i > 3. At large Nc , this expression hits the unitarity bound A = 0 around 
Qym = 0 very quickly. In other words, the operators with the anomalous dimension fl5.8p cannot 
be studied in the ’t Hooft limit. 

The relation between the divergence of the lowest eigenvalue in (15.81) and the pathological 
behavior of YY operators found in |Hl [9] is unclear. The problem in the naive planar limit 
has disappeared by including all non-planar corrections (i.e. by identifying the ground state of 
integrable spin chain with boundaries as a holomorphic-antiholomorphic operator). However, 
there is another operator whose dimension diverges in the large Nc limit with all non-planar 
corrections included. 

Let us recall that there are other situations where the integrability methods such as the gen¬ 
eralized Liischer formula or the mirror TEA disagree with perturbative field-theory calculation. 
Not all the disagreements are related to non-planar effects. One of the most puzzling examples 
is the single-particle state tr {XZ) in the /3-deformed theory [7T1172] . The held theory calcula¬ 
tion shows that this operator is protected owing to the so-called prewrapping effect. However, 

L EE 2 (mod 4), consider C>hoi with = tr {YY) - i tr {X,Xj)- 
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the TBA computation based on a naive asymptotic Bethe Ansatz neglecting the prewrapping 
effect predicts a divergent answer. It is not known how to incorporate the prewrapping effect 
in the integrability methods. 

5.2 Level-crossing and level-pairing 

Below we regard one-loop dimensions as locally real analytic functions of , instead of col¬ 
lection of real numbers evaluated at integer values of Nc . It allows us to keep track of each 
eigenstate from large Nc to small Nc . 

5.2.1 Level-crossing 

We ask the question when two operators have exactly the same one-loop dimensions at finite 
Nc. 

Consider the characteristic polynomial of the non-planar mixing matrix Mjj . This polyno¬ 
mial factorizes into 

<P(7) = det (M/j - 7(5/j) = ( 7 ), (5.9) 

a 

where is a prime polynomial over C at generic values of Nc . 

The roots of different prime polynomials are unrelated, and nothing prevents the level¬ 
crossing. For example, at L = 6 there exists an apparently non-BPS operator which is protected 
at any values of Nc 0 At L = 8 we find two non-BPS operators whose anomalous dimensions 
are simple functions of Nc . These energy levels do cross with the other energy levels}^ 

In contrast, the roots of the same prime polynomial rarely collide unless Nc is small, which 
can be roughly explained as follows. The roots of the same prime polynomial come from a sub¬ 
matrix whose off-diagonal components are non-zero, and the non-zero off-diagonal components 
keep the eigenvalues separated. 

Let us argue more accurately by writing down the condition that the one-loop dimensions of 
two operators coincide at finite Nc . We denote the dilatation operator at Nc = N, and N, — 6N 
by D, and D, — 6D, respectively. We want to solve the eigenvalue equation at Nc = N, — 6N 
perturbatively around 6D = 0. The operator mixing equation can be written as 

(2). - (5D)V'= 7^^ V'= (2).-7)"^ (5.10) 

If we take the limit 6D —)■ 0, only the states with (2), — y)?/) = 0{6D) contribute to the 
right-hand side. In particular, if the energy levels of two states are sufficiently close, the 
other energy levels are neglected. We denote the one-loop dimensions of the two states at 
Nc = N,, N, — SN by 7,(±), 7(±), respectively. The condition that these two states become 
degenerate at Nc = N, — 6N is given by 

7(+) - 7(_) = \/(e Wi - 122)^ + 4 V 12 V 21 = 0, 7,(+) - 7,(_) = e, (5.11) 

is known that this operator remains protected at two-loop order m- 
^'^We use the words “energy level” and “one-loop anomalous dimension” interchangeably. 
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where SD ■ Oj = Vjj Oj. If V/j is Hermitian, then the level-crossing is possible only if the 
off-diagonal components V 12 = V^*i vanish. The off-diagonal components vanish, for example, 
when two states belong to different irreducible representations of the symmetry group of the 
model. This non-crossing rule is a famous statement by von Neumann and Wigner |50| . 

The matrix V/j is not Hermitian in the non-planar mixing problem of A/” = 4 SYM. Two 
roots of the same prime polynomial can collide and become a pair of complex conjugate roots. 
Of course, any gauge-invariant operators of = 4 SYM must have real dimensions, at least 
to all orders of perturbation in qym ■ It suggests that whenever two eigenvalues collide and 
become complex, the corresponding eigenvectors should be nullified by the finite reduction. 
As a corollary, the roots of a prime polynomial will never collide for > lH 

A few remarks are in order. Since the levels repel each other, the energy levels are quite 
dense for At^ > L 3 > 1 , and the one-loop dimensions can change little as we vary as long as 
Nc > L. Also, the l/N^ expansion of the one-loop dimensions fails to converge at the points 
when complex roots show up. 

5.2.2 Specifying branches 

Let us consider the non-planar eigenvalue problem from a mathematical point of view. 

The dimension of a physical operator is a root of the characteristic polynomial fl5.9p evalu¬ 
ated at integer points of Nc . Since the polynomial is an analytic function of Nc , the dimension 
can be analytically continued at any Nf ., which we call an eigenvalue curve. By keeping track 
of eigenvalue curves, we can see which finite eigenstate is connected to which of the large 
Nc eigenstate. 

Not all large Nc eigenstates can be extended to small Nc . As Nc decreases, a pair of 
adjacent energy levels collide and create a pair of complex conjugate energy levels. In other 
words, if we start from a small Nc theory and increase Nc , a new pair of states are created at 
finite Nc . There also exists a pair-annihilation, where two adjacent energy levels collide as Nc 
increases. The combination of pair-creations and pair-annihilations makes it quite non-trivial 
to keep track of eigenvalue curves as a function of Nc . Some eigenvalue curves draw an S-shape 
trajectory on the ( 7 , Nc) plane via a pair-creation and pair-annihilation, as shown in Figure]^ 
Inside the S-shaped region, an eigenvalue curve is no longer a single-valued function of Nc . 

Mathematically, specifying an operator diagonalizing two-point functions of N" = 4 SYM, 
is equivalent to specifying one branch of the algebraic curve defined by the characteristic poly¬ 
nomial (15.911 . We can define a branch of the algebraic curve algebraically or geometrically. 
Algebraically, a branch is a root of the characteristic polynomial (Nc), where 7 ^"^(A'"c) is a 

single-valued function of Nc. Geometrically, a branch is a connected component of the eigen¬ 
value curves on the ( 7 , Nc) plane among the collection of eigenvalue curves. Although the two 
definitions are closely related, it turns out that an eigenvalue curve is not always a single-valued 

^®This fact can also be used as a consistency check on the computation of the non-planar mixing matrix. 
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Figure 6 : The S-shape trajectory of eigenvalue curves based on L = 8 data. 


function of H In short, the geometric definition is more useful than the algebraic one at 
finite N^. Then, how many connected components are there? How precisely can we specify an 
operator of SU{N^) A/" = 4 SYM? 

To answer this question, we must define the meaning of “connected”. It is convenient to 
exclude a point or points where almost all branches meet. In our definition of the one-loop 
dilatation operator fll.ip . we should exclude the points ( 7 , iVc) = (± 00 , 0 ) and study the curve 
in the region > 0. Then, each branch can be specified by prescribing their large N^. behavior, 

^( 7 ,Ar,) = 0, hm y(iV,)=7j“) + iVA'7i“^ + iVA^^“^ + ... (Y, > 0). (5.12) 

Nc^+00 

We define a connected curve by identifying the adjacent branches which collide at finite W. 

Alternatively, we can also use the rescaled eigenvalue 7 = 7 /W- Then, the large N^. 
dimensions can be read off from the asymptotic slope of the curves in ( 7 , W) plane, and we 
exclude the points ( 7 , W) = (oo,oo)|^ More importantly, the 7 ’s stay finite around W = 0, 
and each eigenvalue curve can be smoothly extended to W < 0 . 

The reader may be upset because one cannot define gauge-invariant operators for W < 1 • 
This difficulty can be circumvented by replacing U{N^ gauge group with -|- k\k). The 

one-loop mixing matrix does not change if we modify tr(l) = to str(l) = For a 

sufficiently large fc, the supergroup theory is not subject to any finite W constraints and one 
gets non-unitary AdS/CFT correspondence realized by ghost D-branes [731 [ZH. In this setup, 
the dilatation eigenvalues at negative W are well-defined, and they can be complex due to the 
loss of unitarity. 

Now let us have a closer look at the rescaled one-loop dimensions at negative W. 

^®Moreover, this definition depends on the choice of coordinates on the (7,77) plane. 

^^The large W zero modes asymptotes to (7,77) = ( 0 ,oo). We should also exclude these points. 

^®We neglect the decoupling of U{1) for simplicity. The A? = 4 SYM scalar $ with SU{Nc) gauge group 
should be extended to the scalar $ with PSU{Nc + k\k) gauge group by imposing tr ($) = str ($) = 0 . 
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If we take the limit Nc —)■ — oo, we should recover the spectrum of planar dilatation oper¬ 
ator, neglecting the flipped sign. The invariance under Nc —)■ —Nc can be understood as the 
invariance under the interchange of M <H- in the U{M\N) gauge group. Another explanation 
is that the’t Hooft limit of SU{Nc) theory is universal in the sense that one cannot discern if 
Nc is positive or negative]^ 

Let us give yet another proof that the dilatation spectra at Nc = ±oo are identical. From 
the block structure of the dilatation operator, one can show that the characteristic polynomial 
is compatible with the following Z 2 symmetry: 

^( 7 , AT,) =^( 7 ,-TV,) 77 ^( 7 iV„iV,)=^((- 7 )(-iV,),-iV,). (5.13) 

To explain the block structure, let us take monomial multi-trace operators as a basis of the 
Hilbert space and collect the operators together according to their trace structure. The matrix 
elements of the dilatation operator within the same trace structure (i.e. block-diagonal parts) 
are of order A)?, while those between different trace structures (i.e. block-off-diagonal parts) 
are of order A:""^. Now recalling the definition of the characteristic polynomial, 

^( 7 , TV,) = ^ sign(cr) (TVfi„(i) - 7 5 i^(i))(TlT 2 ^( 2 ) - 7 ^ 2 ^( 2 )) • • • {Mdaid} - 7 <^rf( 7 (d)), (5.14) 

(reSd 

one finds that, for each a, all terms come with an even power of TV,. Hence the identity (I5.13p 
follows. Note that this identity does not imply that all eigenvalues are expanded in powers of 
TV“^, because there may be a pair of eigenvalues obeying ■j± = a± b/Nc + c/Nc -f-.... 

In terms of the rescaled eigenvalues, the above identity induces an automorphism 


I : (7,Wc) >-7 (-7,-/Vc), 


(5.15) 


which shows that the eigenvalue curves and their mirror-images are identical on the (7, TV,) 
plane. In general, the automorphism t relate different eigenvalue curves. To see this, consider 
the relation between the highest and the lowest eigenvalues. 


~ (lowest) = -^(highest) ^ . [^(lowest) _ 


(5.16) 


Our data in Appendix |B.2| shows that the outermost eigenvalues never collide with the adjacent 
eigenvalues, and remain real for real TV, . Thus, the two energy levels are disjoint on the (7, TV,) 
plane. Incidentally, the highest energy state in the 50 ( 6 ) sector is the product of Konishi 
operators, whose one-loop dimension at large TV, behaves as 7 = 3L -|- TV,“^ L{L — 2)/S + .... 

Generally, by using t we can group together the low-energy and high-energy eigenvalue 
curves as 


^(n-th lowest) = n-th highest) j'lyCfi-th lowest) 


(5.17) 


Thus, the upper (or lower) half of the spectrum is redundant. We can throw away a half of 
the eigenvalue curves, or a half of the operators in f/(TV, + k\k) SYM theory by using l. We 
call it level-pairing. The level-pairing structure is evident from Figure where we plotted the 
rescaled eigenvalues for TV, G M at L = 4, 6. The plots for L = 8 is shown in Appendix |B. 2.3 


^®The sign is important to distinguish the SO{Nc) and Sp{Nc) W = 4 SYM theory [7511751177] . 
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Figure 7; The rescaled one-loop anomalous dimensions of so (6) singlets of length L = 4 (upper) 
and L = 6 (lower) at hnite N^. Left hgures show the eigenvalues on the real plane Imp" = 0. 
Right hgures are an aerial view showing both Re 7 and Im 7 . We included the real eigenvalues 
which are not continuously connected to the large W data on the real plane Im( 7 ) = 0 in both 
hgures. The eigenvalues do not collide oh the real plane because of the dimensionality; two 
straight lines in more than two-dimensional space usually passes over or under the other. 
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Here we list some properties of the level-pairing. First, l maps a 1/iVc-protected states to 
itself, because 7 = '^N^. is a straight line on the ( 7 , iVc) plane. Second, to hnd the level-pairing 
one must consider all multi-trace operators, even in the planar limit. Third, the level-pairing 
preserves so ( 6 ) charges. 

The automorphism l induces a non-trivial map between the conformal dimensions. It maps 
A = L -|- |^Ac 7 -|- to i o A = L — |^Ac 7 -|- (9((7ym)- Then one needs to analytically 

continue Nc to —N^. along the same eigenvalue curve to obtain A' = a.c. (to A) = L + -\- 

0{gYM)- If is not easy to quantify the relation between A and A' explicitly. 

Let us make a few speculative comment about the automorphism l. The Z 2 symmetry (I5.13p 
should persist at higher orders in gyu because this is a corollary of the large Nc factorization 
(see |78| for a review). Then we can ask whether the level-pairing can be interpreted as certain 
duality on AdSs x superstring in the planar limit. Such a transformation should act on 
multi-string states on AdSs x S^, and map a low-energy state to a high energy state. It should 
also satisfy other expected properties, namely to map a low-energy state to a high energy state, 
and a 1/Ac-protected state to itself with the same so (6) charges. 

Since t is a relation between multi-string states, a worldsheet duality (such as the worldsheet 
T-duality [7S1 EHl EH IH2]) cannot be the precise counterpart of t. Nevertheless, it may happen 
that a collection of multi-string states can be described by a single-string coherent state if the 
number of strings is sufficiently large. For example, if the automorphism l corresponds to the 
worldsheet T-duality, then the product of Konishi operator in the limit L —)■ cx) would be dual 
to a hoop-like string]^ 


6 Discussion and Outlook 

In this paper, we studied the spectrum of the non-planar one-loop dilatation operator among 
50 ( 6 ) singlets by computing the matrix elements of operator mixing explicitly up to L = 10 . 

At large Nc we considered the submixing problem, which concerns how the large Nc degener¬ 
acy is lifted by 1/Nc corrections, and observed interesting patterns. First, the 1/Nc corrections 
to the dimensions of the large Nc zero modes start appearing at the order 1/A/. Second, the 
coefficient 72 is always non-positive. Third, the operators with a different number of traces do 
not mix at large A^, and those with the same number of traces do mix. Fourth, the submixing 
density satishes the projective commutation relations. 

We have given some expressions of 72 in terms of correlation functions and derived the 
submixing Hamiltonian. Using the correlator expression of 72 , we have shown that 72 is non¬ 
positive in the large Nf approximation. We studied in detail L/2-trace operators and the 

infinitely winding hoop string potentially suffers from two types of instability: a' corrections and gg 
corrections. The leading ol correction to the energy of a hoop string is complex [53], which represents the 
instability of a hoop shrinking into the north or south pole. It is likely that the correction is also complex. 
Since the hoop intersects with itself infinitely many times, each piece of the hoop can recombine itself to string 
bits without consuming energy. The resulting state may be dual to the Konishi product. 
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double-trace operator, which is the simplest operator to work because it does not submix the 
other large zero modes. We also estimated how 72 depends on the operator length L. 

From the AdS/CFT point of view, the negative sign of 72 can be interpreted as follows. 
Multi-trace operators can be regarded as multi-string states on AdSs x S®, and the product 
of BPS operators in A/" = 4 SYM correspond to the product of BPS string states. In general 
a product of BPS states is not protected by supersymmetry. When multi-string states start 
interacting, they attract each other and form a bound state due to gravitational interaction. 
The negative sign should also be related to the causality constraint on the AdSs x S® side [7]. 
To identify the precise relation, we need to clarify two issues. First, we studied so(6) singlets 
which are different from higher-spin operators of [ 3 . Second, we found operators which carry 
small positive anomalous dimensions at 0{1/N^) or higher. 

When the rank of the gauge group is comparable to the operator length, we are studying 
determinant-like operators with all non-planar corrections taken into account. In AdS/CFT, the 
YY double-determinant operators should correspond to a pair of giant- and anti-giant-graviton 
D-branes. Such DD-hrane configuration is unstable and should decay into the vacuum of 
AdSs X S®. At weak coupling, we find two interesting operators of L ~ 2Nc whose behaviors are 
similar to those of the DD-hrane system. The holomorphic-antiholomorphic operator can be 
regarded as non-planar completion of the YY double-determinant operator whose dimension is 
protected. The lowest-energy eigenstate in the so (6) sector has a negative anomalous dimension, 
which can be regarded as the non-planar ground state. We conjecture that the anomalous 
dimension of the lowest-energy state diverges to —00 in the double scaling limit p] 

Furthermore, we have clarified the structure of level-crossing. When adjacent energy levels 
collide, the corresponding eigenstates should become null to prevent complex eigenvalues. We 
found it quite non-trivial to keep track of the operator dimensions as we vary N^. As a by¬ 
product, we discussed there is a natural pairing between different energy levels, particularly in 
the non-unit ary U{M\N) theories. 

Our methods are mostly based on brute-force computation, whose power is limited only 
to relatively small L. The dimension of the mixing matrix grows factorially with respect to 
L, and we encounter a challenging problem even numerically. It takes a long computational 
time to study the property of a huge matrix, and the results become less reliable due to the 
accumulation of numerical errors. 

It is therefore important to look for another approach to the finite Nc problem. In Section 
we related the submixing to the four-point functions of (products of) BPS operators, which will 
be an interesting future direction. Examining the integrability in higher-point functions can be 
insightful, and rewriting the correlators in the Mellin space may also be useful 15^ 157] . 

Another promising approach to work on the finite Nc physics is to exploit group represen¬ 
tation theory. Operators can be labeled by a set of Young diagrams, and the operator mixing 

is not clear what the gravity dual of the lowest-energy state with L ~ 2A7 is. The dual object might be 
a composite of another D-brane system, e.g. |84]. RS thanks Nadav Drukker for discussion on this point. 
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problem is neatly described by group theoretic quantities. Yet, the study of the operator mixing 
problem in the 50 (6) sector has not been very active so far because of the complexity caused 
by the flavor structure. The 5o(6) singlet sector is much simpler than the full 5o(6) sector. As 
we have seen in this paper, the singlet sector contains interesting physics that have not been 
observed in the 5u(2) sector. It would be one of the next directions to look closely into the 
50(6) sector at hnite Nc- 
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A Notation 

Let us dehne real scalars of W = 4 SYM by 

dim G dim G 

(«>.)• = , (•*.); = E 

A=1 A=1 

where G is the gauge group, a = 1, 2,..., 6 and i, j = ^,2,..., Nc. As for G = U{Nc), the 
symbol <F differentiates $ as 

(k)i (tsA = (A.2) 

This rule can be used to compute the one-loop mixing matrix of SU{Nc) theory as well by 
imposing the traceless condition tr (<Fa) = 0 when constructing a basis of operators. The 
splitting-and-joining rules follow straightforwardly: 

I-1 I-1 

tr (A<l>a)tr (<F6R) = tr (A<FaS$fe) = 5^6 U (A) tr (R). (A.3) 

The dilatation operator and conformal dimensions can be expanded in series of 

Stotal = Stree + ®one-loop + ■ ■ ■ , A = Ap + Ai + . . . , Ai = 7- (A.4) 

^^The dimension of Konishi multiplet is A = Ag -I- 6 ^ ^ -I- .... 
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The non-planar dilatation operators in the scalar sector at one-loop is 

2)one-ioop = ^ : (^-^tr[<h^,<h„][cE)”^,<|)”] - ^tr (A.5) 

which commute with all so (6) generators. It also commutes with the projection operator to a 
union of eigenspaces 


[-Px, Stotal] — 0, Px — 1 — It/’/) (V’/l 5 Stotal IV’/) — ^ IV’/) • 

I&H 


The first term of flA.5p acts on multi-trace operators as 


tr : tr ($,A) tr (<h,S) = tr 

tr : tr {^aA^^B) = tr (A) tr (5<hb<h,) + tr (5) tr (A<h,<hb), 

tr [<I>^<I>„<|)"cE.-]; tr ($,A) tr (<hfeS) = tr {BA^a^, + 

tr : tr (<h,A<hb5) = tr (A) tr (5<h„<h,) + tr {B) tr (A<hb<I>,), 


and the second term as 


(A. 6 ) 


(A.7) 

(A. 8 ) 

(A.9) 

(A.IO) 


tr : tr (<h,A) tr (<h,S) 

tr [<h^$”<h^<h"]:tr(<h,A$fe5) 
tr : tr (<h„A) tr {^^B) 

tr [<hA<f>”<f>"]:tr(<h„A$fe5) 


2 5„6tr(<I)^A$^5), 

2 5ab tr (<I)mA)tr 

5ab tr [^m^m{AB + BA)], 

6ab {tr (A) tr ($^<1)^5) -h 


(A.ll) 

(A. 12 ) 

(A.13) 

(5)tr(<hA7l)}. (A.14) 


A.l Notation for operator mixing 


To highlight the structure of the non-planar operator mixing, we introduce a formal parameter 
Nf so that the global symmetry becomes so{Nf). The A/" = 4 SYM corresponds to Nf = 6 . 
There can be several ways to define the dilatation operator at general Nf , but the difference is 
not much important at the fully non-planar level. We use the simplest generalization; namely 
the operator identical to fll.ll) except that the flavor indices run m,n = 1 , 2 ,..., Nf 

The matrix of operator mixing is defined as follows. We use monomial multi-trace operators 
as the basis, and denote them by {Oj}. By applying the one-loop dilatation fll.ip to them, we 
obtain the mixing matrix 

'D-Oi = MijOj. (A. 15) 


Then we look for the eigenvector of the form V’o = Cai Oj , which implies Mjj Caj = 7 Cai ■ Note 
that this mixing matrix is the transpose of flA.lSp . and in Mathematica the eigenvectors are 


Another possible generalization is 


r»T 


1 


■^one-loop 

which is integrable at large TV, 


4 tr [$.,$„] [ 4 .-, 


tr [$^,4>"][$„,4>" 
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given by Eigenvectors [TransposeQM] .Table [0 [i] ,{i ,d}]. The matrix elements of operator 
submixing will be defined in the same way, namely 

H:^-Oi = {M,^)ijOj. (A.16) 

Explicit matrix elements are computed via (I3.30p . 


B Foundations of finite Nc calculation 


B.l Finite Nc reduction 


When the rank of the gauge group is smaller than the operator length, gauge-invariant operators 
become linearly dependent. This is a well-known phenomenon, and the linear relations are 
called finite W constraints,^ Any finite Nc constraints are written as the identity for an 
antisymmetric tensor 


0 ^[n*2--'*jVc+i] ’ 1 ) 2 ,..., Nc)■ 


(B.l) 


As Nc decreases, the dimension of the Hilbert space of states shrinks by finite Nc reduction. 

The two-point functions at tree level can be diagonalized at any Nc by using group-theoretical 
bases uni na [la HZ]. In these bases, finite Nc constraints are typically written as ci{R) < Nc, 
where Ci {R) is the number of rows of a Young diagram R. A set of finite Nc constraints in the 
so{Nf) scalar sector was given in [17] . In general there are several ways to express the finite 
Nc constraints in a given sector. (This is the same as that there are several orthogonal bases 
in a sector.) The reason can be clarified from the existence of some sets of conserved charges 
at tree level [63] . 

At the loop level, there is no freedom in choosing the basis in which the two-point func¬ 
tions are orthonormal, except for degenerate cases. Still, the finite Nc constraints reduce the 
dimension of the Hilbert space of gauge-invariant operators. When the finite Nc constraints 
are imposed, we should find either (i) the eigenvalues remain unchanged, or (ii) the eigenvector 
becomes null. Written explicitly, this means 


25one-loop ”0/ 'yi'4^1 


finite Nc 


® one-loop "07 0/ 0 


(B.2) 


where is the eigenstate subject to the finite Nc constraints. 

Let us make a few remarks. First, if an eigenvector ipj becomes null, all correlators involving 
ijji also vanish, Second, if an eigenvector becomes null at Nc = h, then it remains null for 


Nc < h — 1, as follows from fIB.ip . 


A simple way to determine which eigenvector survives at finite Nc is to substitute Nc x Nc 
real traceless matrices to the fields . In other words, we regard multi-trace operators as 

^^This is a non-perturbative effect on the string theory side, and called the stringy exclusion principle m- 
^®See [91] for an example. 
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GL(A'"c)-invariant polynomials of the matrix elements (*ha)p • When W decreases by one, we 
remove the last row and column of the matrix elements. Schematically, this can be depicted as 
removing gray region of the following matrix; 


$ 


a — 



Nc 

(<ha)ll = -5^(<hak 

i =2 


(B.3) 


By exploiting this idea, one can associate a Young diagram to each eigenstate as follows. 
Instead of removing the last row and column, we can rescale the entries in the last column and 
the last row by ejsf^. The diagonal element (<I>a)are rescaled by . Then we measure how 
an eigenstate scales in the limit; -0 ~ t 0. We can continue this procedure as 


lim ... lim lim 

ejVc— 


-knc-1 

^Nc-1 


^-K,Nc 

-Nc 


^ = 0{el- 


), 


(B.4) 


until h = 2. This assigns to each eigenstate -0 a sequence of integers (^ 2 , ks, ..., kl) which is a 
partition of L. 

The hnite-iVc constraints put restrictions on the matrix elements of non-planar operator 
mixing. Let Vh be the space of the GL{Nc) -invariant polynomials of ($a)p which vanish for 
Nc < h. In other words, h is the smallest integer satisfying Kj = 0 for all i > h. Each Vh 
is closed under the dilatation operator. However, the mixing matrices are generally not block 
diagonal on the basis of {V/i}. To see it, expand the eigenvector which becomes null at h as 


q(W) v^ + J2 Vj , (vi e Vh , V, e V 2 \ Vh). (B.5) 

* j 

The null condition 'ip{h) = 0 gives Cj{h) = 0. Since the coefficients are non-trivial functions of 
Nc , they do not identically vanish. 

That said, the same reasoning allows us to remove some elements of the mixing matrix. If 
we use monomial multi-traces as the basis of operators, then the elements of mixing matrix are 
at most linear in N~^] see flA.7p - flA.14p . This property remains unchanged unless we rotate the 
basis by a A'^c -dependent matrix. Then, the condition Cj{Nc = h) = Cj{Nc = h — 1) = 0 implies 
Cj = 0, which excludes off-diagonal elements from Vh' to Vh for h' < h — 2. 

Here is a side remark. In Liouville theory, the one-point function on a torus is equal to the 
four-point function on a sphere with degenerate helds 1^ IM] . It would be interesting to 
seek for similar relations for the correlation functions of the null states in A/" = 4 SYM. 


B.2 Spectral data 

We summarize the basic properties of the hnite Nc spectrum, and the spectral data of one- 
loop anomalous dimensions for the scalar so(6) singlet operators at hnite W with length L = 
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4,6,8,10. These results are obtained by explicit computation using Mathematica. The data 
include plots of the hnite spectrum, the submixing matrix and the basis of operators. The 
eigenvalues of submixing matrices are given at general values of Nf for L < 8. 


B.2.1 L = 4 


There are four so{Nf) singlets at length L 
general Nf are given bj' 


4. The matrix elements and operator basis for 



( ^ 

-2 

2Nf 

Nc 

— \ 

Nc 


( tr 

Mij = 

-2 

2-2Nf 

Nc 

Nf + 3 

2Nf-2 

Nc 

2 

Nc 

0 

Nf+l 

Nc 

2 

, Oj = 

(02i02l 022^^2) 

tr(0ii0z2)tr(0i2^n) 


\ Nc 

12 

Nc 

0 

2Nf) 


\tr(0ij0ijtr(0i20i2)/ 


(B.6) 


which is consistent with [HI] at Nf = 6. The spectrum is shown in Figure The eigenvalue 
curve is drawn by keeping track of the large Nc eigenvalues down to small Nc as long as 
they remain real-valued. The complex eigenvalues are not shown, and some of them have 
null eigenvectors. Whenever the complex eigenvalues show up, the corresponding eigenvectors 
become null |95) . 

There is one large W zero mode, so the submixing matrix is one-dimensional: 


- = ( 20 ), 0° = CijCij , 


(B.7) 


where Cij is the symmetric traceless tensor dehned in 02 .201) . The l/W corrections to this 
operator for general Nf are given by 


(Nf + 2)(N, - 1){N, - 3) 
N,N2 


0(N-‘) ~ -20JV-^ 


(Nf = 6 ). 


(B.8) 


a computer-friendly format, Mu is 

{{4.-2,-2*Nf/Nc,2/Nc},{-2.3+Nf,-2/Nc.Nc-(-l)+Nf/Nc},{2/Nc-2*Nf/Nc,-2/Nc+2*Nf/Nc,0,2},{-12/Nc,12/Nc,0,2*Nf}} 
One can derive submixing matrix elements and various finite Nc plots from this data as presented below. 
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Figure 8: The dimensions of 5o(6) singlets of length L = 4 at hnite Nc . The upper hgure shows 
the one-loop anomalous dimensions, and the eigenvalues corresponding to null eigenvectors are 
denoted by x. The lower hgure shows the sum of tree and one-loop dimensions as a function 
of Nc and the’t Hooft coupling A = Ncgyu- 
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B.2.2 L = 6 


There are fifteen 5o{Nf) singlets at length L = 6. The mixing matrix and operator basis for 
general Nj are given by 

{{6,3,-6,0,0,(3*Nf)/Nc,3/Nc,-6/Nc,0,6/Nc-(3*Nf)/Nc,3/Nc,-6/Nc,0,0,0}, 

{0,4+Nf,0,-4,2,2/Nc,2/Nc+(2*Nf)/Nc,-6/Nc,0,-4/Nc,3/Nc-Nf/Nc,2/Nc,0,0,0}, 

{-2,1,6,-2,0,-6/Nc+Nf/Nc,Nc"(-1),4/Nc-(2*Nf)/Nc,2/Nc,-6/Nc,2/Nc+Nf/Nc,4/Nc-Nf/Nc,0,0,0}, 
{0,-3/2,-2,6+Nf/2,0,(-2*Nf)/Nc,3/(2*Nc),-2/Nc,2/Nc+Nf/(2*Nc),2/Nc,Nc~(-1)+Nf/Nc,-4/Nc,0,0,0}, 

{0,3/2,0,-3,3+(3*Nf)/2,-6/Nc,15/(2*Nc)+(3*Nf)/(2*Nc),0,-3/Nc,0,0,0,0,0,0}, 

{2/Nc,2/Nc+Nf/Nc,-4/Nc,(-2*Nf)/Nc,Nf/Nc,3+Nf/2,3/2,-2,0,0,0,0,-2/Nc,Nc"(-l)+Nf/(2*Nc),l/(2*Nc)}, 
{0,12/Nc,0,-24/Nc,12/Nc,0,3+2*Nf,0,-2,0,0,0,0,-2/Nc,Nc~(-1)+Nf/Nc}, 

Mij = {-4/Nc,-2/Nc,6/Nc-(2*Nf)/Nc,-2/Nc+(2*Nf)/Nc,2/Nc,-2,0,4,l,0,0,0,-(Nf/Nc),Nc~(-l),0}, 

{0,0,-24/Nc,24/Nc,0,0,-2,0,4+Nf,0,0,0,0,(-2*Nf)/Nc,2/Nc}, 
{6/Nc-(3*Nf)/Nc,6/Nc,-12/Nc,6/Nc+(3*Nf)/Nc,-6/Nc,0,0,0,0,6,3,-6,0,0,0}, 

{-4/Nc, -2/Nc-Nf/Nc, -4/Nc, 4/Nc, 6/Nc+Nf/Nc, 0,0,0,0,0,2+Nf, 0,0,0,0}, 

{-6/Nc,(3*Nf)/Nc,12/Nc-(3*Nf)/Nc,-6/Nc,0,0,0,0,0,-6,3,6,0,0,0}, 
{0,0,0,0,0,-6/Nc+(3*Nf)/Nc,3/Nc,6/Nc-(3*Nf)/Nc,-3/Nc,0,0,0,0,3,0}, 

{0,0,0,0,0,24/Nc,-2/Nc+(2*Nf)/Nc,-24/Nc,2/Nc-(2*Nf)/Nc,0,0,0,0,Nf,2}, 
{0,0,0,0,0,0,36/Nc,0,-36/Nc,0,0,0,0,0,3*Nf}} 


and 

{ri/} — ■^1+^0^^^(^22 013 022^012^13)! 1^(021021023 022 0^2 013)! ll" (021 022 021 023 022 023 ) ! 1^(021 021 022 023 022 023)! 

(0*1 0*1 0*2 0*2 0*3 0*3 ) ’ (0*1 0*2 ) (0*3 0*1 0*2 0*3 ) ’ (0*3 0*3 ) (0*2 0*1 0*1 0*2 ) ’ (0*1 0*3 ) (0*2 0*1 0*2 0*3 ) ’ 

(0*3 0*3 ) (0*1 0*2 0*1 0*2 ) ’ (0*1 0*2 0*3 ) (0*3 0*1 0*2 ) ’ (0*3 0*1 0*i ) (0*2 0*2 0*3 ) ’ (0*1 0*2 0*3 ) (0*3 0*2 0*1 ) > 

tr (0i^0i2 )tr( 0 i 3 0 ijtr( 0 i 2 0 i 3 ), tr( 0 i 2 0 i 2 )tr( 0 i^ 0 i 3 )tr( 0 i 3 0 ij, tr( 0 i^ 0 ijtr( 0 i 2 0 * 2 )tr( 0i3 0i3 )}• (B.9) 

The spectrum is shown in Figure]^ This mixing matrix is consistent with j62], which can be 
shown by computing the characteristic polynomial ^( 7 ) = det (M/j — y^/j) and substituting 
7 —)■ 27, iVc —)■ l/i/, A^/ —)■ 6. 

There are two large Nc zero modes. The submixing matrix for general Nf is given by 

3(Af/-l-4) [3Nj+32Nj-60Nj-132Nf-S) 

4(V/+2) (v2+12Af/-|-14) 

0 

(B.IO) 

The double- and triple-trace zero modes do not submix. At Nf = 6, the double-trace mode 
has 72 = —3675/61 and the triple-trace mode has precisely zero dimension for any Nc |62| . 

We also studied operators other than so (6) singlets. No negative modes (72 < 0) are found 
in the su(2) sector at L = 4, and the lowest eigenvalue in the so(6) sector at L = 4,6 is same 
as the lowest eigenvalue of the so(6) singlets. 



0 

)(^/-2)(jV/+4) 

4:N f 


o° = 


‘ ^ij^jk^ki , 
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B.2.3 


L = 8 


There are 71 so(iVj) singlets at length L = 8. The mixing matrix and operator basis for general 
Nf are given by 



and 




■Ai 3 'Ai 3 ^i 4 ). 

, tv iPi^Pi^Pi^Pi^P.^Pi^Pi^Pi^), 

<t>i2<t>i3<t>i4)' 


<Ai 4 ^i 3 .Ai 4 ). 

, tviPi^Pi^Pi^Pi^P^^Pi^Pi^Pi^), 

<Ai 3 <Ai 3 <Ai 4 ) 

, tr (Pi2 <Ai4 ) tr (pi^ pi^ pi^ <Ai2 <Ai3 <Ai 4 



, tr ((Ai2 Pi^ ) tr (<Ai3 <Ai4 ^14 <Ai2 '#'i3 </>i4 

<l>i2<Pi3) 

, tr (cAig cAij ) tr (<Ai2 '/’i4 'Aii <Ai2 >/'i3 'Ai4 

<^i2<P^2) 

, tr ((Ai4 (Ai4 ) tr (pi^ pi^ pi^ pi^ pi^ pi^ 

<Ai 3 <Ai 4 ) 

, tr (cAij Pi4 Pi2 ) tr (pi^ pi^ <Ai2 Pi^ 4‘i4 

<At 3 0 , 4 ) 

, tr ((Ai2 Pi2 4>i2 <l>i4 ) tr (pi^ pi^ Pi^Pi^ 


, tr (cAij Pi2 Piv Pi2 ) tr (pi^ pi^ Pi^ Pi^ 

to 

, tr {pi^Pi^Pi^Pi^)^^ (Pi^pi2pi2pi2 


tr (4>i^ 4>i^ 4>i2 ) tr (.4>iv 4‘i2 <l>ii <l>i2 ^is)’ (‘t'ii <i’ii '^>i2 ) (‘^•3 ‘^•1 '^•2 '^•3 ‘^•2 ‘^•4 '> ('^*3 ‘^•2 ‘^•3 ( 

tr (<l>i2 <l>ii 4 ‘ii </>i2) tr ((p i2 <l> i i <l> i ^ <l> i 4^) ^,3 c#)i3 ), tr (<#)i2 </)i3 ) tr {</>i^</>ij j ), tr ( 

tr (<^i3 <^i3 <Ai2 '#'i 4 ) <^^3 '#'i2 ‘^>2 ('^^4 '^•2 "^^3 )’ ('^*3 '^*3 '^•2 '^^4 '> “^^2 ( 

tr <t>i2<t>ii <l>i2) tr {(Aij'#’14'^>13'^>14 ), tr ) tr {<l>i^<t>i^) tv <t>i2<t>i2)< tr {<^^3 ) tr <^^3 ) tr (</>i2 <#>i4 (Ai^ <Ai2 ) ■ 

tr (cAij 0i2 ) *t ('#'i4'#’t3) tr ■Ai2‘^»3)' *t W*1 '^^4 ) tr (i/'i4'#>.2 ) ('^is'/'ii <Ai2'^*3)’ tr (lAiiiAia) tr ((pi^tpi^) tr {4)i^<l>i^ i/'i2>/'i3)> 

tr (<Ai3 <Ai4 ) tr {<Ai^ <Ai3 ) tr ((Ai^ <Ai2 </>ii 'Ai2 )’ tr (<Ai^ <Ai3 ) tr (,Ai4 'Ai2 ) ('^^4 <^^2 '^^3 ("^n "^^3 ) ('Ai2 '^il ’^i2 "^^3 )’ 

tr (<Ai3<Ai3) tr {rl>i^<t>i^) tv (<Ai2 <Aii <Aii <Ai2 ): tr {<Ai3 (Aig ) tr (<Ai4 iAi4 ) tr (cAij <Ai2 ■#>12 ) ’ *t (’^i3’#'i2 ) ('^n <^^2 ) ('^M ’^i4 ) ■ 

tr (<Ai3 0i2 ) tr {<Ai3 cAij <Ai4 ) tr {<Ai^ (Ai^ <Ai2 )■ tr {■Ai4 <Ai4 ) tr (lAi^ .Ai2 '#>t3 ) tr (iAi3 4>i2 </>ii ). tr ((Ai3 (Ai^ ) tr (pi^ ) tr {<pi^ (pi^ pi^), 
tr (<Ai4cAi4) tr {<Aij (Aij (Aij ) tr (rpi^rpii 4 >i 2 )t tr {iAi3 <Ai2 ) tr (lAi^ iAi3'Ai4 ) tr ((Ai 4 lAii <Ai2 ). tr (iAi4<Ai4) tr {<pi2pi2 <Ai3) tr ( 0 i 3 lAii <Aii ) , 
tr (<Aii <Ai4 ) tv (<Ai2 <Ai2 ) tr (<Ai3 <Ai3 ) tr ( 0 ^^ , tr (cp^^ ) tr {(Ai^ pi^ ) tr (pi^ pi^ ) tr (<Ai^ pi^ ), tr (pi^ pi^ ) tr (p^^ p^^ ) tv (pi^ p^^ ) tr {pi^ pi 

tr iPi^ Pi4 ) tr ipi^ Pi^)tv (pi^ pi^ ) tr (pi^ pi^), tr {cAij <Aii ) tr (pi^ pi^ ) tr (pi^ pi^ ) tr (pi^ 0i4 )} ■ 


1 ). 

(B.ll) 


The spectrum 


is shown in Figure 
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Figure 10; The one-loop anomalous dimensions 7 and the rescaled dimensions 7 = iVc 7 of 
50 ( 6 ) singlets of length L = 8 at hnite . In the upper hgure, eigenvalue curves do not cross 
for A'";; > T = 8 and continuously connected to the large N^. spectrum. Some of the eigenvalue 
curves have a loose end at small Nc, due to creations and annihilations of complex eigenvalues 
as shown in the lower hgure. 
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Also, there is a pair of exceptional operators at L = 8, which take the form 

^ = tr (0jj0j2) tr [ 0430^2043 , 0 j^ 0 j^]) + ^ tr (043 0^3 042 013 0*4 [ 0 * 3,042 014]); 

1 

18 


^ ( 5 iV, ± V 25 iV 2 + 504 ) , 


and their one-lop dimensions are exactly given by 


7(^^ = ^(33± 


^/ 25 N^ + 504 ' 


Nr. 


(B. 12 ) 


(B. 13 ) 


There are hve large Nr zero modes. The submixing matrix for general Nf is given by 


{ ^sm') ij 


/ (Nf + l)(Nf+6)xii 

3(Nf+2){Nf+4)^10 

2(iVj+l)(Ar/+6)jt2i 

3 (]Vj+ 2 )(JVj,+ 4)4<20 

(jV/-2)(iV/+l)(iV/+6)>.22 

2(jV3'+2)2>ir20 

0 0 


Y 0 0 


0 

2Nj{Nf + l){Nj+6)>t22 
3(iV3+2)(]Vy:+4)j€20 
[Nf + l)(^Nf+6)^23 
2iVj(iVj+2)2,(20 


0 

0 


, , , , 6(N'i-Nf+6) 

{Nf - 6) {Nf + 3) ^ ^ 


4 {Nf - 6) 






(B.14) 


Kii = 20823Af)^ + 1675098Ar)° + 55911035Af^ + 1008163871Af^ + 10662939982Vj + 66380603808V^ + 219710331899Af^ 
+ 178394761318Af| - 1267103763092V^ - 4642599661272Ar| - 6369708861792Af/ - 3243010514688, 
xio = 1971Af^ + 158898Af| + 5341569Vj + 97989373Nj + 1077096339Af| + 7325782233V^ + 30674558032V| 

+ 75803848404V| + 99816658464V/ + 53440944192, 

K21 = 124V| + 2271Nf + 8447V| - 15528V^ - 88236V/ - 81648, 

K22 = 63V| + 1153V| + 4522V| - 2208V/ - 12960, 

K23 = 163V^ + 2459V^ + 876Nf - 63460V| - 94656V| + 136800V/ + 217728, 

K20 = 37N^ + 727N] + 3714V/ + 4536 (B.15) 


These matrix elements are dehned on the basij^ 

{f^ijklC'ijkli Cijkl^ij^kly ^ijk^ijl^kly \f^ij ^jk ^kl ^li\ ; ^ij ^ji ^kl ^Ik^ ; 

where [CijCjkCkiCn]' is the following combination of two quadruple traces, 

[C4J Cjk Cki C;*] tr (043 042) tr (042 043) tr (0*3 044) tr (0*4 0*3) 


(B. 16 ) 

(B. 17 ) 


+ 


3 

Wf 


tr (0*3 044) tr (0*4 041) tr (04^043) tr (0*3 0*2) 


— tr ( 0*4 0 * 4 ) tr ( 0*3 0*2) tr ( 0*3 0 *i) tr (0*2 0 * 3 ). 


By diagonalizing the submixing matrix at iVj = 6, one hnds 


{ 7 ^~ {- 101 . 276 ,- 66 . 6009 ,- 44 ,- 6 . 89905 , 0 }, (B. 18 ) 

where 72"*^ are 72 for m-trace operators. The operator corresponding to 72^^"*" is not protected, 
and receives corrections at 0 {Nr^) with a positive sign. 


^^The ofT-diagonal elements of Mgm depend on the normalization of C. Here we use the convention: = 

li 'Ea(^Se tr (^*,,(1) ■ • ■ ®i<,(o) “ (^avor contractions). 
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Figure 11; Spectrum of the one-loop dimensions among the 5o(6) singlets of length L = 10 at 
hnite Nc. Each eigenvalue curve is continuously connected to the large Nc spectrum. 


Degeneracy of non-zero modes. There are six degenerate positive eigenvalues at L = 8, 
namely {6,7,8,9.5,12,14}. The eigenvalues receive hnite Nc corrections at the order 1/Nc, 
and the large Nc expansion of these modes can submix with the operators having different 
number of traces. They make a contrast to the properties of the large Nc zero modes flB.ldp . 
For example, the eigenstates corresponding to the eigenvalue 6 take the form 


~ 05,3 -|- 03 , 3,2 , O 2 ~ 04 , 2,2 , O 3 ~ 02 , 2 , 2,2 , O 4 ~ 05,3 -|- 03 , 3,2 , 
where Oli,L 2 ,... E a multi-trace operator with length {Li,L 2 ,... }. 


(B.19) 


B.2.4 L = 10 

There are 469 so{Nf) singlets at length L = 10. The mixing matrix and operator basis for gen¬ 
eral Nf are shown in the attached hies AncillaryNegative.nb and AncillaryL10Data.txt. 
The spectrum is shown in Figure [TT|p^ To highlight the structure of the spectrum lying in the 
middle, we plot the eigenvalue density in Figure [T^ 

There are 11 large Nc zero modes, as shown in flD.lOD . There are no hnite-A^c zero modes. At 
Nf = 6, the submixing matrix has two zero modes corresponding to combinations of quadruple- 

^®We have not checked potential level-crossing with exceptional eigenvalues due to the complexity. 
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Figure 12: The eigenvalue density of the mixing matrix length L = 8 (left) and L = 10 
(right) at hnite . The vertical axis shows the number of real-valued eigenvalues in each bin 
[7 — A 7 / 2,7 -|- A 7 / 2 ). The bin width A 7 is adjusted for each L. 


and quintuple-trace operators. The one-loop dimensions of two operators are 0{N^ and 
positive. The submixing matrix at Nf = 6 is given by 


( ^sm) ij — 


132.315 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 0 0 
87.4095 36.7996 40.8718 
34.5075 50.803 0 


25.2692 

0 

0 

0 

0 

0 

0 

0 


77.1909 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

36.6545 

-2.64168 

-12.0246 

-14.1498 

-10.4861 

0 

0 


0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

-0.0000105607 -12.6725 -21.2528 0 

1.97708 2.62507 2.22822 -9.00914 

-6.72677 59.5308 -34.3732 30.6522 

-1.17112 -23.5379 40.8351 5.33674 

-13.2197 7.32801 6.00623 60.2393 

0 0 0 0 

0 0 0 0 


26 


0 

0 

0 

0 

0 

0 

0 

0 

0 


16.5923 0/ 


which is block-diagonal. These matrix elements are dehned on the basis 

0° = {C^, 3 X 5 X 2 X C^}, 


(B.20) 


(B.21) 


where p x C"* means p linear combinations of the large Nc zero modes with m-traces, and 
explicitly given in flD.lOp . The eigenvalues of the submixing matrix is 


72 = -{132.315, 125.673, 83.7144, 67.7793, 66.8384, 42.3268, 26, 22.892, 5.41638, 0, 0}. 

(B.22) 

Only numerical values are shown in (lB.2np because their precise forms are quite complicated. 
For example, the eigenvalue for the double-trace mode is 


72 = 


36834342860563635266164905898354615349151994033997 

278383379192667725097781339828443592271084673096 


-132.315. 


(B.23) 
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Figure 13; The one-loop dimensions of the 5o{Nf) singlet large Nc zero modes with L = 8,10 
at Nc = 10. Blue, orange, green curves correspond to Nf = 4,6,8, respectively. 


B.2.5 L = 12 

There are 4477 so{Nf) singlets at length L = 12. We have not made a detailed analysis of the 
spectrum due to a huge amount of computation involved. There are 34 large W zero modes, 

0° = {C\ 5 X C^, 14 X 10 X 4 X C^}, (B.24) 

in the notation of flB.2ip . 


B.3 Spectrum at general Nf 

We consider the spectrum of submixing matrix for 5o{Nf) singlets at general Nf, which will 
highlight the special properties of A/” = 4 SYM from the non-planar spectrum of other gauge 
theories. The value Nf = 6 corresponds to A/" = 4 SYM. 

Let us outline two properties of the spectrum of the submixing Hamiltonian at general Nf . 
First, 72 decrease as Nf increases. 


d'f2 

dNf 


< 0 , 


(72 <0, Nf ^ 6). 


(B.25) 


It shows that the leading 1/Nc corrections to the large Nc zero modes can be positive when 


Nf < 6 as shown in Figure 13 


Second, all eigenvalues of the submixing matrix scale as {Nf/Nc^ in the large Nf limit 


72 , 

7 = ^ + 0(iV7) ~ c(L) 


(72 < 0), 


(B.26) 
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The limit Nf ^ 1 does not commute with W ^ !> so we consider W ^ Nf ^ 1. 
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Figure 14: Plot of 72,00 = 72 (6/-^/)^ for large zero modes at Nf = 100. The dashed lines 
represent the Nf = 6 data. 


with some coefficient c{L). This scaling behavior can be seen from the submixing matrix at 
L = 4, 6, 8 given in (jR 
not bad. Figure 


14 


fIB.lOl) . (IB. 141) . Moreover, the large Nf approximation is numerically 
shows that the rescaled eigenvalues 72,00 = 72 (6/A^/)^ do not deviate a lot 


from the Nf = Q data. 


C Correlation functions at large Nf 

We explain the details of the computation in Section 


C.l A proof of {'iJjq'DiiIjo) < 0 at large Nf 


In this subsection, we will show that the first term of 72, {'ipo'^No), is always non-positive by 
making use of the limit iVy 3> T ^ 


Since SqV’o = 0, we can consider ('^o^one-ioopV’o)- For convenience, we rename each term of 
the one-loop dilatation operator by 


2)one-loop = ^ : (C.l) 

where 

dW = -itr[<F^,<l>„][<|)™, <!>"], ^ _ltr [<F^,cE)T[<I>^,cE)”], (C.2) 

and 

Jji2) ^ y^(21) ^ y^{22) _ _ + ^tr (C.3) 

In general k’o is given by a linear combination of the large W zero modes as k'o = ■ Because 

off-diagonal expectation values (* j) cannot have the leading power with respect to both l/N^ 

and 1/iV/) consider only diagonal expectation values ('0o*^®i'0o*^)- 
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Suppose that acts on 


(p = tr tr (<h(bi<hb 2 ■ ■ ■ ■ ■ ■ <hc,)), (C.4) 

where one derivative in acts on one of the <h’s in the first trace and the other derivative 
acts on one of the <h’s in the second trace. We denote the remaining traces of V’o by (p, i.e. 
-00 = W- 

We can find that the leading behavior of iVj in the correlator {'ipo'DI'lpo) is . The 

leading power arises if the following two conditions are satished; (a) the two derivatives in 
2)one-loop act on two ‘hb/s, giving = Nj and replacing them with <hm’s, (b) when a $6^. (or 
<hm) is Wick-contracted with a $ 0 ^,, the other ^bj (or <hm) is also Wick-contracted with the other 
, giving Because the factor cannot arise from the action of only the 

action of may be considered in what follows. At large Nf, the traceless part of '00 does 
not give the leading contribution. 

When the conditions (a) and (b) are satisfied, the correlator can be reduced as 

(w(2)(^)0) = ((p(2)(p))((p(p) ^ ((^(2)(^)). (C.5) 


At the last step, we performed Wick-contractions in (00) and simplified the flavor indices 
using the h’s. In the reduced correlator, all flavor indices are contracted. Then, the dilatation 
operator makes a single trace from (p, removing two ‘hb/s and putting two <hm’s. More precisely, 
put two <hm’s away from each other while put them next to each other, as shown in 
(A.11) and (A.13). When two <hm’s are adjacent in the reduced correlator, we have a smaller 
symmetry factor than the non-adjacent case. This difference in symmetry factor is the reason 
the correlator (0o2)i0o) is non-positive. 

The following is an example. Consider a reduced correlator 


-(: tr($,,d), 3 ---d),Jtr($,^<l>, 2 ---<l>,J : . tr (<l>b,<I>b 2 ■ ■ ■ <hbj tr ($b,<I>b 2 ■ ■ ■ <hbj :)• (C.6) 

s 

This is a case of p = 0 and q = 0. It has the following two contributions 

J(; tr ($„, ■ ■ ■ <I)„J tr ($,, ■ ■ ■ <I>,J : ; tr (^b, ■ ■ ■ <I>bJ tr (<l>b, ■ ■ ■ <I>bJ :) 

= -Nf{-. tr ($,, ■ ■ ■ $„Jtr ($„, ■ ■ ■ <I)„J :: tr (^^^b^ ■ ■ ■ 'hb.^'m'hb^ ■ ■ ■ <hbj :) 

= -2sN^^-^Nf^+ 0{N^^-^), (C.7) 

The symmetry factor 2 s in can be explained as follows. When ai in the first trace is Wick-contracted 
with the first m, oi in the second trace must be Wick-contracted with the second m. To keep the planarity, we 
should Wick-contract the first 02 with the first 62, and the second 02 with the second 62, and so on. We say this 
that (oi, 02 • • • , Og) are Wick-contracted with (m, b2, - ■ • ,bs). Likewise (oi, 02 • • • ,as) can be Wick-contracted 

with (bk+i, ■ ■ ■ ,bs,m,b2, ■ ■ ■ , bk) for any k while keeping the planarity. In this way, we get the factor s. The 

factor 2 comes from exchanging two single traces. 
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ancE3 

-(: tr • ■ ■ $aj tr • ■ ■ <^>aj : : tr ■ ■ ■ <l>bj tr ■ ■ ■ ^bj ■) 

s 

-/V/(- tr ($Q|J^$£12 ■ ■ ■ “^as) tr ($ai*^a2 ' ' ' *^03) • • tr *^^2 ' ' ' *^6s) •) 

= 2iVf-iiV^+i + C>(iVf-2), (C.8) 


where we have kept only the leading term with respect to Nf and Nc- 
correlator (C.6) is negative. 


Because s > 1, the 


C .2 Nf and L dependence of 72 

In this subsection, we will consider the 1/Nf expansion to study the L-dependence of 72 for the 


operators presented in Section |3.3[ General properties of 72 will be first discussed, and then 
concrete examples are shown. 

Recall the expression of 72; 


72 = 


Ni'Do'ijji) 


No'lpo) 


(V’oV’o) 


(C.9) 


As we discussed in the last subsection, the first term behaves as 0{Nf), and as we will see 
below the second term has a more dominant behavior 0{Nj). 

The relation between 'ipQ and "01 is given by the equation Do'ipi + SiV’o ~ 0 obtained in 


(3.21), which can rewrite the second term as 


NiDNo) 


(C.IO) 


{i>oi>o) Noi>o) 

Focusing only on the leading of the 1/Nf expansion, we find that SiV’o can be expanded as 


NfFi + NfGi + Oil), 


(C.ll) 


where 0(1) is multi-traces without the factor Nf. Fi is a multi-trace that does not contain two 
adjacent matrices with the flavor indices contracted, and Gi is a multi-trace that does contain 
such matrices. They appeared in the image of and respectively; see the second lines 


of flC.7l) and flC.8l) as well as (C.15) for the double-trace operator. The 0(1) comes from the 
action of Od). 

The planar dilatation operator has the form S)o = I — P + G/2, where /, P and O are the 
identity, the transposition and the contraction acting on nearest neighbor matrices. Only the 


^^The symmetry factor (C.81 comes from the Wick-contractions between (01,02 
only. It can be more manifest if we write the correlator as 
(: tr ($ai$a2 • - to,) tr (^a^ • • • $a3$ai) :: tr ($„$&, • • • 


s) and ( 


ro, 02,- • • 


N) 


^) :)• 


rj/Nc in ( 3.211 is not important at the leading order of the 1 /W expansion in 72, so we can use DoV'i + 
TiNq ~ 0. 
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contraction can bring a factor of Nf. In order to determine the form of ipi, we need to solve 
Do'ipi ~ —Nf{Fi + Gi). Since C ■ Fi Gi and G ■ Gi NfGi, the solution is roughly written as 

~-aiV/F, -/3G, + ■ ■ • , (C.12) 


where a and /? are A^/-independent coefficients. Substituting these equations into the form of 
72 , we have 


72 = -N^a 


{F.F,) 

(V’oV'o) 


+ 0{Nf), 


(C.13) 


where the ratio of the two-point functions is independent oi Nj. The A^| behavior of 72 is 
consistent with (B.26). It is emphasized that Gi in (C.ll) does not matter for the evaluation 
of 72 at the leading order. 


Our next goal is to estimate the L-dependence of the correlators. The double-trace operator 
and L/ 2 -trace operators are studied below. 

Consider the singlet double trace operator (3.36). We will use the property of the symmetric 
representation {a{a)\P[L/ 2 \\b) = {a\aP\^L/ 2 ]\b) = {a\P[L/ 2 ]\b)■, and 


{c\P[L/2]\d) ~ YJ2 ^ {cL/2-l\P[L/2-l]\dL/2-l) 

~ j^Nf{cL/2-l\P[L/2-l]\dL/2-l), (C.14) 

where cl/ 2 -i = (ci, ■ ■ ■ ,cl/ 2 -i), and the factor ^ comes from the normalization of the pro¬ 
jector. 

Acting with on ^jJo, we have 
= -^tr 

~ ~ (4’m4’ci ■ ■ ■ ^CL/2-i^m^di ' ' ' 

= -■^^/(cL/2-l|-P[L/2-l]|^4/2-l)tl' (4’m*hci ' ' ' ' ' ' ^dL/ 2 - 1 ) 

= ■ ■ ■ <hc,/2_7), 

= -\N,h, 

= ^tr 

L 

~ ($m‘hm‘h(ci ■ ■ ■ 4’ciy2-i)^(ci ■ ■ ■ ^Ci^/ 2 - 2 )) 

= (C.15) 
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The multi-traces fi and /2 belong to Fi and Gi respectively. Note that the flavor indices are in 
the irreducible representation [L /2 — 1 ] ® [ 1 ], where [p] represents the rank-p symmetric traceless 
representation. 

We have to solve the equation Do'ipi -f ~ 0. Suppose that the solution is given by 




(C.16) 


where a is an A^j-independent positive constant that may depend on L. Other terms are 
denoted by ■ • •, expecting that they do not change the L-dependence at the leading order of 
the 1/Nf expansion. We will leave it as a future problem to determine ipi. We then have 

(/i/i) 


72 


-a 


{i’di’o) 


+ 0[N,). 


The normalization was computed in (3.38), which is evaluated at large Nf as 
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(V’oV’o) = y DimL /2 


N 


Lj2 


2 (L/2)! " ’ 


and we also have 


(/i/i) 


N\ 


L/2 


(L/2 - 1)! 


nL 


to obtain 


72 




(C.17) 


(C.19) 


(C.20) 


(C.21) 


Therefore the double-trace operator scale as 72 ~ aLp- 


The next study is about L/2-trace operators such as CijCjkCki ■ ■ ■ Cmi and CijCjiCki ■ ■ ■ Cmi- 
These multi-traces have the same color structure, but the flavor indices are contracted differ¬ 


ently. As we mentioned in Section 3.3, such multi-traces can be classified by the partitions of 
L/2, and we write them by where T expresses a partition of L/2. 

In general, a submixing eigenstate "i/o is a linear combination of tpQ ’’s. This situation 

(1~'\ ('Y'\ 

simplifies at large Nf , where the two-point functions are orthogonal {ip^ Vo ) oc Vt' and the 

iT') 

submixing Hamiltonian acts diagonally on . Then 72 is expressed by 


(T) 

72 = 


(v-rso-i/r’) 


+ 0(N,). 


(C.22) 


The dimension of the rank-L/2 traceless symmetric representation of so{Nf) is 

1 (L/2+ W-3)1 


Dimi/2 = 


{Nf-iy. (L/2)! 


L(y _l)(L + y-2) 




(L/2)!’ 


(C.18) 


where we have approximated using Stirling’s formula for Nf ;g> L. The result (C.191 can also be obtained by 
considering only Wick-contractions that satisfy the condition (b) in the last subsection. 
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For = Cif.uCki ■ ■ ■ 
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we obtain 


^{ 2 i)^(L/ 2 ) ^ _^iv^tr (<F^<l>,<F^<Ffc) CkiCis ■■■Cu = 

The two-point fnnctions can be compnted 

{44^44"^) ~ ln4^n4 

(^(L/ 2 )^(L/ 2 )) _ 2(1 + 5l,^)N4^N^ 

to obtain 

^(i/2) = - J(1 + &,4)0 (i/2 ,L + 0(Nf), 


(C.23) 


where ol(i^i 2 ) is determined by the equation 

2 ) o^/’i+ 2 ) 1-00 ~ 0 ^ 


,(L/ 2 ) T AL!"^) , 

01 - ^ a{Li2)Lfl '' + ■•■ 


(C.24) 


(C.25) 


(C.26) 


The two expressions (C.21) and ( C.25| ) agree at L = 4 as expected. The equation flC.25p shows 
that the L/2-trace scale like 72 ~ a(^L/ 2 )L as a function of L. 

Consider more general L/2-trace operators Fp/ 2 ) _ ^[Li/ 2 ]^[L 2 / 2 ] ... (j[Lp/ 2 ]^ where 

L = and we assume Li 7 ^ Lj {i 7 ^ j) for simplicity. We then have 


jj{21)^{Li/ 2,--- ,Lp/2) _ 


_ (j[ni/2]/j[i2/2] . . . q[Ls-i/2\ 

9 / 


X (^tr 


■ ■ ■ C 

0304 '^O.Ls/20‘1 


^ l^Fs+l/2] _ _ _ ^[Lp/2] 




(Ll/2,-,Lp/2;s) 


(C.27) 


and 


(0(Li/2,-,Lp/2)^(Li/2,-,Lp/2)) _ ^^/V+S^yL^ 

^^(Li/2,-,Lp/2;.)^(Li/2,-,Lp/2;s)^ _ 2(1 + 5^^ 0+, . . . 

= 2(1+4.,4)0+''"A'/. 


(C.28) 


where Lt = + 1+2 ■ ■ ■ Fp. Substituting these two-point functions into the expression of 72 , we 
will obtain 

r2 t 2 ^,,_ . / 4 ^-iT;s) r{T-,s)\ 


72 


(T) 


+ E,i 7 (T,.)(/!'" 7 r-'> 


{+'v+> 


+ 0(iVi) 


Af; 


~~4 + ^Ls,4)a(r;s) + 0{Nf), 


(C.29) 


where tt('r;s) is another positive constant in i/j^ = (-^// 2 ) Xls'^('L; s)-^s/i^’*^ + ■■■; which is 
determined by the equation 2)o0i + 2)i0o = 0. We leave it for a future problem to determine 
the precise L dependence. 

Introduced the notation CA\ = CaiaadoaosCajai • ■ -Capai- 
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D Relation to Mandelstam Variables 


We pursue the coincidence between the number of so(6) singlet large W zero-modes in Table 


in Section [27^ and the number of the completely symmetric polynomial of Mandelstam variables 
under certain conditions in Table 2 of [00] • This fact allows us to construct an explicit basis 
for the latter. 

We defined as the number of so (6) singlet large W zero modes with length L. These 
zero-modes can be written explicitly in terms of the traceless symmetric single-trace operators 




{}, 

{Cij Cij }, 

{^ijk ^ijk 1 ^ij ^jk ^ki} •> 

{^ijkl ^ijkl 1 ^ijkl ^ij ^kl •> ^ijk ^ijl ^kl •> ^ij ^ji ^kl ^Ik i ^ij ^jk ^kl •> 


(D,l) 


corresponding to {Z 2 , ^ 4 , Zq, Zq} = {0,1, 2, 5}. 

Another series of positive integers is computed in [601, which is the number of the completely 
symmetric polynomial of Mandelstam variables of degree K, subject to massless momentum 
conservation, for a sufficiently large number of particles n. If we denote this series by Z '^^, they 
are given by {Z[, Z 2 , Z'^, Z'^} = {0,1, 2, 5}. The Mandelstam variables are defined by 


«« = {Pi +P^f = 2pi -Pj, (i,j = 1,2, 


(D.2) 


and they satisfy 


^ij ^ji 1 


sp= 0 , (=°- 

j=l \ i=l / i<j 


(D.3) 


A basis of the completely symmetric polynomials of {sp} can be given as 




E 


■S(T(l)(T(2)'Sa-(2)(T(3)'5a'(3)(T(l) 



(1),7{2) 1 


a 






a a 

■®t7(l)(T(2)'5(7(3)(T(4) 5 So-(i)o-(2)'So-{2)o-(3)So-{3)(t(4)'So-(4)(t(1) |, 

a (7 


(D,4) 


where we sum a over the permutation group of n-th order Sn ■ 
imposed after the summation. 


The constraints flD.3p must be 
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The relation between flD.lD and flD.4p can be explained graphically. Note that we neglect 
dimensionality constraints; namely the hnite Nf constraints in (ID.ip and the Gram determinant 
constraints in (ElP 


We begin by the 5o(6) singlet large W zero modes. Since is symmetric traceless, 

a flavor index of should be paired with the flavor index of another Since is 

symmetric, the position of the flavor index inside is irrelevant. When we contract m indices 
of and we draw m lines in between as 


C ■ C ■ = 



Similarly, the zero modes with length six are expressed as 


C ■ C ■ = 


C C- ■ c ■ = 


A 


and with length eight as. 


A' A 


I- - I 


I - I 


(D.5) 


(D.6) 



(D.7) 


Let us turn to the completely symmetric polynomials of Mandelstam variables Given a graph 
representing the large W zero modes, we label each “single-trace” by i = 1, 2,..., n. Then for 
each line connecting the Gth and j-th trace, we associate Mandelstam variable Sij , as 




Sl2 


Sl2 


= s ^2 + permutations. 



(D.8) 


which gives the first line of flD.4p . Similarly, a basis of the completely symmetric polynomials 
of degree three is expressed as 


= 5^2 + permutations. 


A 


= -Si2 S23 S31 -|- permutations. (D- 9 ) 


^®The latter comes from the linear relations among {pf} when n is greater than the spacetime dimensions. 
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One can check that the length eight graphs flD.7p reproduce the polynomials of degree four in 

dElD. 

This pattern continues. At L = 10 we have Z'^ = 11. The correspondence can be seen from 


and 


^ijklm ^ijklvn i ^ijklm ^ijk ^Ira i ^ijkl ^ijkm ^Im i ^ijkl ^ijm ^klm i 

^ijkl Cij Ckm Clm i C-ijk ^ijl C-km ^Im i C-ijk C-ilm ^jl C-km j f'pTc f'p'fc Om Om 5 

f'p'fe ^Imk ^ij ^Im ) f'p ^ji ^kl ^Im ^mk 1 ^ij ^jk ^kl ^Im ^mi ) 


(D.IO) 


S12 + • • ■ , ■sL'Sia + • • • 5 S^2'^13'523 + • • • , Si2'Si3'®23 + • ■ ■ , 

'512S13S14S34 + ■ ■ • , S12S13S24S34 + • • • , S12S13S14S23S24 + • • • , 'Sl2'^34 + ■ ■ ■ ; 

'512S13S34 + • • • , Si2S34S35'S45 + • • • , S 12S23S34S45S51 + . . . , (D.ll) 


where +... means the sum over the permutations Sn ■ Note that the all terms add up with 
the same sign. 


It is not obvious to prove the linear independence when the massless momentum conserva¬ 
tions fID.Op are imposed. For example, there are three candidates of the symmetric polynomials 
at degree two, 


-^2,1 — S 12 + ■ ■ • 5 -^2,2 — S 12 S 13 + ■ • • 5 P2,3 — 512534 -|- . . . . (D.12) 

It turns out that P 2,2 = —P 2 ,i and P 2,3 = 2 P 2 ,i when the constraints fID.Op are imposed. This 
can be shown, for example, by applying the condition Ylj Sij = 0 repeatedly, until there are 
no indices which appear only once. We checked that the polynomials written in our “graphical 
basis” are linearly independent up to P < 5. It also follows that our basis is complete thanks 
to Z 2 K = Zj^ . 

We also emphasize that it is non-trivial to construct an explicit basis of the completely 
symmetric polynomials of {s^}. A candidate is given by taking multi-traces of the matrix of 
Mandelstam variables. 


P{k, £,...)= tr {Sk^Sk 2 ...) tr ...), Sk = 


( 0 5^2 
’?2 0 




5 2 n 


sL... 0 


, (A:>0). (D.13) 


/ 


\'5ln ^2n 

These polynomials are manifestly invariant under permutations, but highly redundant even 
before imposing the massless momentum conservation. It is not clear how to take a linearly 
independent set. 

When we turn on masses, the momentum conservation becomes 

2 


j=i \i=i 


Pi I = 


5ij — 2 Mi , 


(D.14) 


i<j 


2=1 
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and the number of independent completely symmetric polynomials increases. For example, 
three degree-two polynomials in flD.12p are no longer proportional. The corresponding gener¬ 
alization in A/" = 4 SYM will be to consider U{Nc) gauge group instead of SU{Nc). 

It is interesting to see how the above relation can be proven for general K, and whether the 
correspondence can be extended to massive cases. 


E Polynomial notation for multi-traces 


A new concise notation for multi-trace operators is invented to reduce the computational work¬ 
load significantly in Mathematica. We call it polynomial notation, which can be used to describe 
any gauge invariant operators of A/" = 4 SYM. 

We associate a polynomial to each multi-trace operator. The basic idea is to describe the 
position inside a trace at which a given SYM field appears, rather than to describe which SYM 
field appears at a given position inside a trace. 

First, consider single-trace operators. If a field F stands at the p-th position inside the 
trace, we write F^, then sum over all fields as 

tT{XYZZY^Z^) o P{X + Y^ + Z^ + Z^ + Y^ + Z^), (E.l) 


where P means that this polynomial is defined modulo cyclic permutation of a trace, 

P{f + g^ + --- + + z^)=P{f + g^ + --- + y^ + z). (E.2) 

For later purposes, we rewrite this equation as 

P = T- P, (E.3) 


where T is the shift operator defined by 


T-P 




CL L 

(a < 0) 


a 

(1 < a < L) . 

(E, 4 ) 

a — L 

(L < a) 



Next, we generalize this notation to the operators with flavor indices contracted, like XiX\ 
for the u(3) singlets and for the so(6) singlets, where Xi = ^ 2 i-i + i ^ 2 i ■ If Xi is at the 
p-th position and Xj is at the g-th position, we write 


tr(...w...x/...) P{<jYf+ ...). 


(E.5) 
This polyno- 


The so (6) singlets can be expressed by imposing the relation ~ ( 
mial notation is significantly simpler than the usual one because no flavor subscripts are used 
anymore. To see this, consider the following u(3) singlet operator of length four. 


tT(x,,x,,xlxl)=tT(x,,x,,xJxl) o 


t vt 






(E,6) 


keep 4> to avoid confusion between and 
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On the left-hand side, we wrote two identical operators generated by the permutation of flavor 
indices t (v(i) 0 ( 2 ) • • • 0 (l/ 2 )) with a G 5^/2- As L increases, the order of the 

permutation group grows factorially. It becomes quite cumbersome to keep removing redundant 
elements generated by In this sense, the polynomial notation is factorially simpler to 

describe single-trace operators than the conventional notation, 
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Second, consider multi-trace operators. If a field F stands at the p-th position inside the 
m-th trace, we write F^, then sum over all fields as 


tT{YX,)tT{YXj) o P{Y, + Y, + <P'i<p;), 


(E.7) 


where (t^(t\i means that we contract Xj at the p-th position in the m-th trace and Xj at the 
g-th position in the n-th trace. We impose the cyclicity condition P as 


P = Tm ■ P for each m. 


(E.8) 


where Tm is the shift operator flE.411 for the m-th trace. 

The polynomial notation is related to a graphical representation of multi-trace operators as 
follows: 





= P U(h+ cYcb + 


(E.9) 





P ( 0102 + 0102 + 02 ^ 2 ) ) 




0^ 


(E.IO) 


where the black dashed lines represent the SU{Nc) color traces, and the blue solid lines represent 
so(6) flavor contractions. The middle equality in the last line comes from the cyclic symmetry 
of the trace. In short, the polynomial notation describes the configuration of the blue lines. 
Note also that a similar graphical method was useful in finding the relation to Mandelstam 
variables in Appendix 

We generated all so(6) singlet multi-trace operators up to L = 12 using Mathematica, and 
implemented the action of the non-planar dilatation operator in the polynomial notation. The 
details are given in the attached file. 
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